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ABSTRACT 


It  is  the  purpose  of  this  study  to  consider  the  post-buckling 
behavior  of  an  axially  compressed  thin  circular  cylindrical  shell.  To 
this  end,  a  simple  model  of  the  post-buckling  configuration  has  been 
devised.  The  model  required  the  consideration  of  three  parameters  and 
minimization  with  respect  to  two  parameters.  The  model  is  shown  to 
produce  results  comparable  to  current,  more  elaborate,  investigations. 
As  with  the  current  solutions  to  the  problem,  the  results  of  the 
analysis  presented  in  this  thesis  are  strictly  qualitative. 
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CHAPTER  1 


INTRODUCTION  AND  DISCUSSION  OF  THE  LITERATURE 

1 o 1  Introduction 

This  thesis  will  consider  the  problem  of  determining  the 
elastic  post-buckling  behavior  of  a  thin  circular  cylindrical  shell 
under  axial  load,  a  problem  first  defined  in  its  present  form  by 
von  KArmAn  and  Tsien  in  1941.  It  can  be  shown  that  the  mathematical 
model  of  von  KArmAn  and  Tsien  is  governed  by  two  nonlinear  partial 
differential  equations,  the  derivation  of  which  is  based  on  a  neglecting 
of  the  nonlinear  reference  to  the  tangential  displacement  components 
in  the  geometric  relations.  The  desirable  procedure  in  effecting  a 
solution  would  be  to  solve  these  differential  equations  subject  to 
the  imposed  boundary  conditions  and  other  constraints,  but  this  has 
not  yet  proven  possible.  Von  KArman  and  Tsien  initiated  a  consideration 
of  the  solution  through  the  principle  of  virtual  work.  Many  investi¬ 
gators  have  since  continued  the  approach  of  von  Karman  and  Tsien,  and 
the  essential  difference  between  their  work  and  that  of  von  KArman  and 
Tsien  has  been  the  addition  of  more  parameters.  This  extension  of  the 
number  of  parameters  from  the  three  von  KarmAn  and  Tsien  used  to 
twenty  or  more  has  been  made  possible  through  the  use  of  modern  electronic 
computers . 
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Nearly  thirty  years  have  passed  since  von  Karman  and  Tsien’s 
basic  work,  yet  it  cannot  be  stated  that  substantial  progress  has  been 
made.  The  post-buckling  load  has  been  lowered  by  successive  investi¬ 
gators,  but  a  less  often  discussed  fact  is  that  the  theoretical  geometrical 
results  deviate  considerably  from  the  experimental  geometry.  The  use 
of  additional  parameters  has  by  now  resulted  in  a  post-buckling  load 
which  is  less  than  the  experimental  result.  This  would  seem  to  indicate 
a  major  theoretical  difficulty  in  the  procedures  used  in  present 
investigations  of  this  problem,  for  one  expects  that  a  prescribed 
continuous  form  of  buckled  shell  would  be  stiffer  than  the  actual 
buckled  shell,  and  should  therefore  result  in  an  upper  bound  to  the 
post-buckling  load.  Until  such  time  as  these  problems,  the  discrepancy 
between  the  theoretical  and  experimental  geometry  and  the  theoretical 
and  experimental  loads,  have  been  resolved,  the  theoretical  results 
cannot  be  said  to  be  a  firm  foundation  for  design;  rather,  the  designer 
must  have  recourse  to  experimental  data. 

Because  of  the  poor  agreement  between  experiment  and  theory, 
it  was  felt  that  a  simpler,  less  involved,  analysis  could  be  found 
which  would  at  least  achieve  results  comparable  to  those  of  the  present 
theory  and  would,  hopefully,  yield  results  closer  to  the  experimental 
values.  To  this  end,  a  simple  model  of  the  post-buckling  configuration 
has  been  devised.  This  model  required  the  consideration  of  three 
parameters  and  minimization  with  respect  to  two  parameters.  The  model 
is  shown  to  produce  results  comparable  to  the  current,  more  elaborate, 
investigations.  As  with  the  existing  solutions  to  the  problem  of 
post-buckling,  the  results  of  the  analysis  given  in  this  thesis  are 
qualitative  and  are  not  suitable  for  design. 
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1 ' ^  Discussion  of  the  Literature 

Before  undertaking  a  discussion  of  the  literature  on  the 
subject  of  post-buckling,  it  is  useful  to  examine  the  experimental 
observations  concerning  the  post-buckling  configuration  of  a  thin 
circular  cylindrical  shell  under  axial  loado  Some  necessary  notations 
are: 


o  =  axial  compressive  stress; 

R/t  =  radius  to  thickness  ratio; 

E  =  Young's  modulus; 

OR/Et  =  nondimensional  axial  stress; 
n  =  number  of  waves  in  the 

circumferential  direction. 

The  experimental  observations  are,  in  part,  as  follows: 

1)  The  cylinder  buckles  from  a  uniformly  compressed  small- 
deflection  configuration  to  a  large-deflection  state  (termed  the 
"post-buckling  configuration") .  The  actual  failing  stress  varies 
greatly,  and  is  probably  dependent  upon  such  factors  as  initial  imper¬ 
fections  and  eccentricity  of  load  application.  Upon  buckling,  the  load 
drops  greatly  to  a  value  termed  the  "post-buckling  load". 


2)  The  large-deflection  configuration  is  composed  of 
diamond-shaped  patterns;  the  diamonds  are  approximately  square,  and 
commonly  occur  as  shown  in  Figs.  1  and  2. 
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Fig„  1  Figc  2a 

Ideal  Buckle  Tier  Buckle 


Helical  Buckle 


Fig.  2c 
Random  Buckle 


The  "Ideal",  shown  in  Fig.  1,  seldom  occurs  in  tests  but  it  serves  as 
the  simplest  mathematical  model. 

3)  Buckling  usually  occurs  over  only  a  part  of  the  shell 
and  is  predominantly  inwards.  Figure  3,  which  is  from  de  Neufville  (1) , 
gives  an  indication  of  the  buckle  size  and  shows  that  n  increases 
as  R/t  increases. 


12 
10 
n  8 


6 

4 

2 

0 


12.8 


J - ! - 1- 

2  4  6 


L  =  cylinder  length 


-I - ! - —i - - — — — 

8  10  12 


R/t  X  10“2 


Figc  3  Variation  of  n  with  R/t 


* 


4)  The  nondimensional  buckling  load  depends  upon  R/t  in 
that  it  decreases  as  R/t  increases „  This  fact  is  illustrated  in 
Figo  4  taken  from  Weingarten,  Morgan,  and  Seide  (2) „  The  wide  scatter 
of  the  experimental  results  is  evident „ 


R/t  X  1CT2 

Fig0  4  Variation  of  Experimental  Buckling  Load  with  R/t 


5)  The  post-buckling  configuration  can  be  made  largely 
elastic  by  a  suitable  combination  of  material  and  R/t  ratio 0 

This  chapter  will  now  consider  in  turn  the  most  pertinent 
theoretical  papers  dealing  with  "large-deflection"  and  "post-buckling" 
theories  of  circular  cylindrical  thin  shells  under  axial  loado  Only 
those  papers  dealing  with  the  static  equilibrium  behavior ,  as  opposed 
to  the  dynamic  behavior  [see  Yao  (3),  Yao  (4) s  Bieniek,  Fan,  and 
Lackman  (5)  ,  Roth  and  Klosner  (6)  ]  s  are  discussed,,  The  papers  reviewed 
are  Donnell  (7,  1934)*  von  Karman  and  Tsien  (8,  1941),  Donnell  and 
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Wan  (9,  1950),  Kempner  (10,  1954),  Almroth  (11,  1963),  Hoff,  Madsen, 
and  Mayers  (12,  1966),  Uemura  (13,  1963),  and  Pogorelov  (14,  1962). 
Reference  will  also  be  made  to  the  works  of  Yoshimura  (15,  1951), 
Jones  (16,  1966),  and  Madsen  and  Hoff  (27,  1965). 

In  reviewing  the  papers,  the  co-ordinates  and  components  of 
displacement  of  a  point  on  the  middle  surface  of  the  shell  will  be  as 
shown  in  Fig.  5. 


Fig.  5  Co-ordinates  and  Components  of  Displacement 

Classical  buckling  theory  [see,  for  instance,  Pogorelov  (14)] 

shows  that 

aR/Et  =  0.60. 

The  fact  that  experiments  failed  to  confirm  this  relationship  led  to 
the  need  to  extend  the  classical  theory.  The  basic  extensions  to  date 
are  those  of  Donnell  (7,  1934),  who  introduced  the  application  of  non¬ 
linear  large-deflection  theory  plus  the  effect  of  initial  imperfections 
and  of  von  Karmdn  and  Tsien  (8,  1941),  who  introduced  the  concept  of 
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a  large-deflection  post-buckling  state.  The  papers  published  since  1941 
largely  differ  from  these  in  the  extent  of  mathematical  computation. 

1.2.1.  Donnell  (7,  1934) 


As  previously  mentioned,  experiments  do  not  confirm  the 


classical  buckling  load;  indeed,  experimental  failure  stresses  range 
from  approximately  0.6  to  0.1  of  the  classical  value.  The  extreme 
variation  of  results  is  indicated  in  Fig.  4.  Donnell  attempted  to 
explain  this  variation  by  the  existence  of  initial  imperfections  which 
increased  rapidly  under  loading.  It  was  assumed  that  the  load-carrying 
capacity  was  reduced  upon  attainment  of  the  yield  stress  at  the  most 
highly  stressed  point.  His  procedure  is  as  follows  and  is  included 
for  completeness,  since  it  contains  the  basic  equations  of  large- 
deflection  shell  theory. 


Following  Donnell,  let 


wi  =  initial  displacement 


W2  =  final  displacement, 


w 


W2  -  wj  =  displacement  due  to  load, 


and  assume  that 


w  “  wj ; 


so  that 


K  +  1  w 
2 


The  unit  strains  in  terms  of  the  displacements  are: 


e 


x 


r  + 
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3v 

3y 


w 

R 


xy 


3u 

3y 


K 

/  3w 

+ 

2 

\3y 

3v  < 

3x 

+ 

K  - 

3w  3w 
3x  3y 


and  the  changes  in  curvature  are: 


x 


and 


32w 

3x2 

32w 

ay2 


32w 


xy  3x3y 

The  nonlinear  terms  in  w  lead  to  the  terminology  "large-deflection". 
The  absence  of  nonlinear  terms  in  u  and  v  should  be  noted. 

The  stress  resultants  are  defined  by: 

Et 


■x 


1-v 

Et 

1-v 

Et 


„  .  e 
2  \  x 


+  ve 


„  .  e  +  ve^ 
2  \  y  x 


xy 


2 (l+v)  xy  5 


M. 


Et- 


x 


My 


12 (1-v  ) 
Et3 

12 (1-v2) 


2s  VKx  +  VKy 


Ky  +  VKX 


and 


M. 


xy 


Et J 

12 (l+v)  Kxy 


and  the  transverse  shears  are  defined  by  the  equilibrium  equations. 

The  macroscopic  equations  of  equilibrium  (i.e.,  the  equili¬ 


brium  equations  in  terms  of  the  macroscopic  stress  tensors)  number 
six,  Donnell  assumed  that  the  equations 


could  be  approximated  by 


3T 

3x 


x 


+ 


9T 

x' 

3y 


=  0 


and 


3T 


3T 

~  -p  ■'  =:  Q 

9x  3y 

These  equations  become  identities  with  the  introduction  of  a  stress 
function,  F$  such  that 
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and 
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By  eliminating  u  and  v  from  the  strain-displacement 
relations,  the  following  constraint  on  e  ,  e  and  e  is  obtained 


xy 


3  2e 


x 
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which,  if  one  notes  that 
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and 


becomes 


t_ 

d 


a2! 

ax2 


a2F 

v  ay2 


xy 


f  V4p 


2(1  +  y ) t  a  2f 
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Donnell  assumed  w  and  used  this  equation  to  determine  F. 


(1. 2.1.1) 


For  w  Donnell  assumed 


t  /  ~  •  2jtx  2ttX  ^  A  _  4ttx 


w  =  -  ^  W  sin  L 


sm 


+  A  cos 


x 


x 


(1.2. 1.2) 


where 

=  axial  wave  length 

and 

L  =  circumferential  wave  length. 

y 

The  first  term  in  Eq.  1.2. 1.2  is  the  classical  solution  and  the  second 

term  is  designed  to  reduce  the  large  deflection  stresses.  This  assumed 

form  of  w  introduces  the  four  unknowns  W,  A,  L  ,  and  L  into  the 

x  y 

problem. 


Returning  to  Eq.  1.2. 1.1,  the  assumed  form  of  w  would 
indicate  that  a  particular  solution  for  F  is  of  the  form 

4tt 


F  =  B  sin  sin  +  C  cos 

x  y  y 


(1.2.1. 3) 


Substitution  of  Eqs.  1.2. 1.2,  and  1.2. 1.3  into  Eq.  1.2. Id  yields 

B  =  B(W,  A)  (1.2.1.4a) 


and 
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C  =  C(W,  A)  0  (1. 2dc4b) 

The  complementary  solution  of  Eq,  1.2. 1,1  is  taken  to  be 


which  corresponds  to  a  uniform  axial  compression. 


The  unknown  A  is  determined  from  the  condition  of  loading 

that 

j>  Tydy  =  0  (1,2,1. 5) 

To  determine  W,  Donnell  employed  the  virtual  work  equation 
3W  =  0  >  (1.2. 1.6) 

where  the  potential  energy,  V,  is  defined  by 

V  =  U  -  f _  T.U.dS  , 

;sT  i  i 

U  being  the  elastic  energy,  S  the  surface  of  prescribed  loading, 

T_^  the  components  of  prescribed  loading,  and  Ik  being  the  displacement 
components  on 


Equation  1.2. 1,6  yields 


P 


(X  +  S)2  X 

X  +  (X  +  S) 

(W  +  Wi  )  KZSZW 

W  8 
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(1.2, 1.7) 


where 
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4tt  2  Rt 
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4n2  Rt 
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1  aR 
P  '  d  Et 

Equation  1.2. 1.7  yields  P  in  terms  of  W,  Wj ,  X,  and  S.  It  should  be 
noted  that  the  quantity  Wj ,  the  amplitude  of  the  initial  imperfection, 
appears  in  this  equation  for  P  because  K  is  a  function  of  W. 
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This  follows  from  the  relation 

W]_  = 


W! 

W 


(K  -  1) 

w  =  - -  w 

2 


Now,  if  a  (i0e. ,  P)  is  prescribed  and  X,  S,  and  are  known  then 
Eq.  1,2. 1.7  defines  W„ 

However,  Donnell  wished  to  include  the  restriction  that 
yielding  occurred  at  critical  points  on  the  cylinder.  He  employed 
the  von  Mises  yield  criterion  for  plane  stress  which  is 
oj  -  a  a  +  a  2  +  3a2  =  Y2 
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this  yield  criterion  becomes 

p  2  =  p2  +  pkxw2  ^  0.25  +  3(2  ~  V-^S  j  +  3XSW2 


X 


(X  +  S) 
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(1.2. 1.8) 


y  Et 

Substitution  of  the  expression  for  P  from  Eq.  1.2. 1.7  into  Eq.  1.2. 1.8 

leads  to  an  equation  for  W,  and  the  critical  value  of  P  may  then  be 

obtained.  Donnell  assigned  L  ,  L  ,  and  W  so  that  his  theory  approxi- 

x  y 

mated  his  experimental  results  (the  procedure  being  described  in 
detail  in  Donnell’s  paper).  Donnell’s  basic  results  are  contained 


in  Fig.  6a,  b,  and  c. 
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Fig,  6a  Theory  Approximation  to 
Test  Results 
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Rig,,  6b  Values  of  X  and  S 
Assumed  in  Theory 


Figc  6c  Values  of  Wj  Needed  to  Duplicate 
Experiments 

Donnell  did  not  show  that 
62V  =  0  , 

which  would  seem  to  be  a  necessary  condition  at  the  buckling  configuration. 
Von  KdrmAn  (17)  noted  that  the  initial  amplitudes  of  Fig.  6c  are  excessive 


' 
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and  would  be  easily  visible,,  Experiments  indicate  that  buckling 
proceeds  without  these  large  initial  deflections.  Von  Kdrmdn  further 
noted  that  in  very  thin  shells,  the  failure  does  not  appear  to  be 
plastic,  since  the  deflections  completely  disappear  when  the  applied 
load  is  removed. 

1 . 2 . 2 o  von  Karman  and  Tsien  (8,  1941) 

In  1939  von  Karman  and  Tsien  (17),  in  a  study  of  the  stability 
of  spherical  shells  loaded  by  external  pressure,  introduced  the  concept 
of  the  post-buckling  load  to  shell  theory.  Considering  a  spherical 
segment  as  shown  in  Fig.  7, 


Fig.  7  Spherical  Shell  Segment 


von  K£rmdn  and  Tsien  found,  by  large-deflection  procedures,  that  load- 
deflection  curves  of  the  nature  indicated  in  Fig.  8  were  possible. 
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Fig  0  8  Illustration  of  Post- 
buckling  Load 


The  envelope  of  the  minimum  pressures  was  taken  to  define  the  post- 
buckling  pressures  under  large  deflections. 

Von  Kdrman  and  Tsien,  applying  Donnell fs  basic  theory  of  1934 
to  cylinders,  found  that  results  similar  to  those  of  Fig.  8  held  for 
cylinders.  There  is,  however,  a  significant  difference  in  the  spirit 
in  which  they  applied  the  equations.  Donnell  assumed  initial  imperfections 
and  attempted  to  find  the  upper  buckling  load  where 
SV  =  0 

and 

62V  =  0  ; 

although,  as  previously  noted,  Donnell  did  not  demonstrate  that  the 
second  equation  was  obeyed.  Von  K^rmdn  and  Tsien,  however,  applied  the 
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theory  (with  K  =  1,  i.er,  no  initial  imperfections)  to  find  the 

post-buckling  configuration  where 
6V  =  0 

and 

62V  >  0 


Experiments  reveal,  as  previously  discussed,  that  the  post- 
buckling  configuration  is  composed  of  large-deflection,  diamond-shaped 
buckles.  This  was  approximated  by  von  KSrman  and  Tsien  as 


where 


m  =  yn 

[For  purposes  of  interpretation  it 

Michielsen  (18,  1948)  and  let 
,  ,  fj. 

f0  +  4  =  §0 


is  more  convenient  to  follow 

f0  =  go  "  f1 

fl  =  2gl 
f 2  =  ^g2  -  81 


which  gives 


w 

R 


=  go 


gl 


cos 


mx 

R 


cos 


R 


g2 


cos 


2mx 

R 


+  cos 


This  is  a  superposition  of 

w0  =  go  i  a  uniform  expansion, 

W1  =  gl  cos  |p  cos  >  the  classical  solution, 
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and 

a  /  tnx  ny  \  /  mx  .  ny  \ 

w2  =  g2i  cos  ^  ^  )  cos  [  T  +  ^  ) 

which  has  nodal  lines  in  a  helical  pattern.] 

Substitution  of  m  =  pn  into  von  Karman  and  Tsien’s  format  yields 


The  quantities  n  and  p  were  prescribed  by  von  Karman  and  Tsien  and 

the  equilibrium  equations  used  were 

j> V  =  3V  3V  = 

3  f  0  3  f !  3f2 

The  equation 


was  used  with  the  constraint  that  v  must  be  periodic,  in  order  to 
show  that  the  average  value  of  T^  was  zero. 


The  nature  of  von  Karman 's  and  Tsien 's  results  are  given  in 
Figs,  9a,  and  9b,  where 
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Fig.  9b  Post-buckling  Curve 


It  is  evident  that  no  minimum  is  realized  for  y  =  0.5;  however,  the 

results  do  indicate  the  possibility  of  post-buckling  states  at  loads 
much  less  than  the  classical  buckling  load. 


Von  Karman  and  Tsien  stated  that  the  jump  from  the  small- 
deflection  state  involves  the  release  of  elastic  energy  and  thus 
explains  the  rapidity  of  the  buckling  process.  The  results  of  this 
thesis,  although  highly  qualitative,  indicate  that  the  elastic  energy 
of  the  large-deflection  state  is  larger  than  that  of  the  small-deflection 
state.  This  statement  is  based  on  the  assumption  that  the  "Ideal" 
shape  of  Fig.  1  is  attained. 
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Figures  9a  and  9b ,  when  applied  to  a  specific  R/t,  indicate 
that  as  n  and  y  vary,  the  post-buckling  load  varies.  The  question 
then  arises  concerning  the  best  representation  of  the  post-buckling 
curve.  One  method  of  representation  might  be  to  take  the  envelope  to 
the  n  =  constant  curves,  derived  according  to  von  Karmen* s  and  Tsien's 
method,  but  there  is  still  the  problem  of  determining  y  0  Legget 
(19,  1946)  suggested  that  for  a  given  y  ,  the  curve  could  be  obtained 
by  equating  the  areas  under  two  successive  n  curves  [see  Fig,  10  from 
Michielsen  (18) ] , 


0  o  5  1.0  1.5 


eR/ 1 

Fig,  10  Calculation  of  Areas 

which,  since 

U  =  Jo(e)de, 

corresponds  to  conservation  of  elastic  energy  between  the  two  states 


Michielsen  (18,  1948)  suggested  that  n  and  y  be  allowed 
to  change  freely  and  so  introduced  the  equations 
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Michielsen  stated  that  his  computations  agreed  with  those  of  Legget 
and  Jones  (20,  1942). 

Thielemann  (21,  1963)  has  noted  that  the  variation  with 
respect  to  n  is  inadmissable  since  non-integral  values  of  n  will 
violate  the  constraint  that  the  displacements  must  be  continuous.  The 
same  objection  logically  holds  for  y,  if  the  length  of  the  cylinder 
and  the  boundary  conditions  are  prescribed .  This  objection  is  usually 
overcome  by  assuming  the  cylinder  to  be  infinitely  long.  It  is  not  at 
all  clear  to  this  author  that  Thielemann’s  objection  is  valid  in  the 
light  of  the  fact  that  the  methods  of  analysis  are  themselves  highly 
approximate.  Experiments  would  seem  to  confirm  that  n  achieves 
preferred  values,  and  one  would  therefore  search  for  a  criterion  to 
determine  these  values.  The  variation  of  n  about  an  integer  could 
be  a  criterion,  and  since  the  discontinuity  induced  is  infinitesimal, 
it  is  probably  a  less  serious  error  than  the  inherent  inaccuracy  of 
the  approximate  method,  and  it  probably  results  in  a  structure  less 
stiff  than  the  real  one.  Michielsen's  (18)  findings  would  tend  to 
confirm  this  idea.  Jones  (16,  1966)  carried  through  Thielemann’s 
suggestion  (i.e.,  to  prescribe  n)  in  detail.  As  the  criterion  to 
determine  n,  he  used  the  integral  value  of  n  which,  for  prescribed  a, 
gave  the  minimum  value  of  potential.  Madson  and  Hoff  (27,  1965)  also 
performed  computations  which  involved  prescribing  n. 

Despite  the  above  objections,  the  variation  of  n  has  been 
used  in  many  papers.  This  thesis  shall  employ  the  variation  of  n 
about  an  integral  value.  Indeed,  in  the  one  example  given  by  Jones  (16), 
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it  appears  that  his  final  equilibrium  value  of  n  =  7  could  have 

been  more  rapidly  computed  by  simply  adding  to  his  equations  the 
equation 

™  =  o 

8n 

and  taking  the  nearest  integer  to  the  solution* 

The  above  points  are  carefully  and  clearly  considered  by 
Mushtari  and  Galimov(39) .  They  point  out  that  the  variation  of  n 
is  permissible  if  the  finite  deflection  extends  over  only  part  of  the 
shell  (as  is  usual  in  shell  experiments) ;  it  is  then  necessary  to 
assume  that  the  part  not  covered  by  the  assumed  finite  deflection 
pattern  (termed  by  Mushtari  and  Galimov  as  the  zone  of  edge  effect) 
is  small,  and  that  its  energy  contribution  is  negligibly  small. 

1 o 2 o 3 o  Donnell  and  Wan  (9,  1950) 

Von  Kdrmdn  and  Tsien  (8)  studied  the  post-buckling  state 
of  initially  perfect  cylinders.  Donnell  and  Wan  extended  these  ideas 
to  cylinders  with  initial  imperfections  and  found  that  the  peak  load 
(not  the  post-buckling  load  which  occurs  at  much  higher  deflections) 
was  reduced  from  the  classical  value.  Donnell  and  Wan  varied  u  and 
n. 


Their  method  is  similar  to  that  of  Donnell  (7),  except  for 
the  following  major  changes: 

i.  The  w  format  is  altered  to  that  of  von  Karman  and 


Tsien  (8) . 


■ 


lie  In  line  with  von  Kerman* s  and  Tsien*s  (17)  criticisms 
of  Donnellys  1934  paper,  the  process  was  assumed  to  be  elastic, 
but  the  computations  were  checked  to  determine  whether  or  not 
the  maximum  stress  had  reached  the  yield  stress  prior  to  the 
attainment  of  the  peak  load. 


with 


Thus  the  form  of  w  assumed  was 

/  ynx  ny  ,  ,  2pnx  „  2ny  ,  ,  ^ 

w  =  at  (  cos  —  cos  — +  b  cos  — - —  +  c  cos  — +  d 

\  R  R  R  R  / 


Wi 


IJL 

a 


w 


o 


To  determine  the  value  of  a-j  ,  Donnell  and  Wan  assumed  that 

a  reasonable  initial  shape  for  a  column  could  be 

clz 

ait  =  — - 


which,  extended  to  the  two  dimensional  surface  of  a  plate,  would  become 

U  1x1y 

ai  =  7Z  ~^T~  ; 

where  lx  =  L^/2  and  1  =  Ly/20  U  is  termed  an  "unevenness  factor" 

and  it2  is  included  to  improve  the  format „  The  form  of  ai  was 
further  altered  to  account  for  the  fact  that  the  process  of  rolling  a 
plate  into  a  cylinder  should  tend  to  decrease  the  imperfection  amplitude 
in  the  y  direction^  hence  a  more  correct  formulation  was  assumed  to  be 


ai 


1-q 


9 


with 
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q  =  0.5  and,  noting  that 


al  =  7T73^ 

The  equilibrium  equations  used  by  Donnell  and  Wan  were 

3a  3b  “  3c  ~  3n  ~  3p  "  0 

The  parameter  d  was  determined  from  the  periodicity  constraint 
on  v.  Donnell  and  Wan  were  able  to  solve  the  equilibrium  equations 
for  the  range 

0*75  <  a/o  .  <  1 

—  el¬ 
and  beyond  this  they  kept  b,  c,  and  p  constant.  The  subscript  cl 

denotes  the  classical  value.  The  major  results  of  Donnell’s  and  Wan’s 

analysis  are  shown  in  Fig.  11. 


Donnell  and  Wan  assumed 


and 


1  =  2ttR/ 2n 
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1  =  1  /p 

x  y 


they  reduced  ai  to 


UR7 


Fig.  11  The  Effect  of  U  on 
Peak  Load 


24- 


As  previously  mentioned,  Donnell  and  Wan  checked  their  computations 

to  see  if  the  maximum  stress  had  exceeded  the  yield  stress  before 

the  applied  load  reached  its  peak  value*  Their  conclusions  were  that 

the  load-carrying  capacity  could  be  reduced,  due  to  elastic  failure 

at  the  peak  or  by  plastic  yielding  before  the  peak  load  is  reached* 

This  more  complete  method  of  calculation  (as  opposed  to  the  1934  paper) 

gave  better  values  of  a^,  as  needed  to  explain  the  experimental 

results.  Using  U  =  0,00015,  it  was  found,  for  instance,  that  at 

R/ t  =  1500,  the  value  of  a^  was  approximately  one.  Thus  these 

results  explain  in  a  reasonable  way  why  a  is  not  usually  reached. 

cl 

They  do  not  explain  the  low  value  of  the  experimental  post-buckling 
load.  The  remaining  papers  to  be  discussed  have  as  their  main  interest 
the  post-buckling  load. 


1  *  2 . 4 *  Kempner  (10,  1954) 

It  was  noted,  in  the  previous  discussions,  that  von  Karman's 
and  Tsien's  1941  work  did  not  produce  acceptable  minimums  in  the  post- 
buckling  range.  They  varied  only  the  three  parameters  d,  a,  and  b  =  c, 
while  n  and  y  were  assigned. 

As  previously  discussed,  Michielsen  (18,  1948)  extended 
von  Karman's  and  Tsien's  method  to  include  minimization  with  respect 
to  n  and  y.  This  procedure  resulted  in  a  nondimens ional  load  of 
0.195  which  occurred  at  a  y  value  of  approximately  0.4.  Legget  and 
Jones  (20)  performed  similar  computations.  Although  Michielsen' s 
result  is  much  less  than  the  classical  value,  it  is  still  considerably 
above  the  experimental  minimum  of  0.06  as  given  by  von  Karman  and  Tsien  (8). 
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Kempner’s  extension  consisted  of  minimizing  the  potential 
energy  with  respect  to  all  five  parameters  a,  b,  c,  n,  and  y ;  he  also 
determined  the  stability  of  the  equilibrium  states  so  found.  The 
procedure  used  to  obtain  the  basic  equations  differs  from  that  of 
Donnell  and  von  Karmanj  therefore,  it  will  be  discussed  briefly. 
Kempner's  geometric  equations,  stress-strain  law,  and  energy  expressions 
were  those  of  Donnell 9  except  that  w  was  taken  positive  radially 
inwards.  The  equilibrium  equations  were  derived  by  variation  of  the 
potential  energy  and  yielded 

32w 


m  3zw  ,  3zw 

DV4w  =  T  -r — i  +  2T  -r  - 
x  ox  xy  dxdy 
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(1.2. 4.1) 


(1.2. 4. 2) 
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3T 
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H - - 

3x  3y 


=  0 


(1.2. 4. 3) 


Equation  1.2. 4.1  is  recognized  as  the  plate  equation,  with  a  correction 
of  1/R  for  the  initial  curvature,  Kempner  enforced  Eqs.  1.2 .4. 2  and 


1.2,4. 3  by  defining 

32F 

3yz 
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and  determined  F  from  the  compatability  equation 
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(1.2. 4. 4) 

Equations  1,2. 4.1  and  1, 2.4.4  are  the  governing  partial  differential 
equations  for  the  von  Kdrman  and  Tsien  mathematical  model.  As  mentioned 
in  the  Introduction  their  analytical  solution  has  not  yet  been  obtained. 
Indeed  because  of  the  neglect  of  nonlinearity  in  the  geometric  equations, 
it  is  probable  that  these  equations  are  not  valid  for  deflections  larger 
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than  t  and  are  therefore  not  applicable  to  the  post-buckling  problem. 
Von  Karman  and  Tsien  (8,  1941)  and  Donnell  and  Wan  (9,  1950)  have  made 
similar  comments  but  this  most  important  fact  seems  to  have  been 
completely  overlooked  by  many  investigators  and  they  continue  to  apply 
the  equations  to  the  post  buckling  state  which  involves  deflections 
of  50t  or  more. 

The  procedure  is  not  sufficiently  general  to  yield  the  correct 
equilibrium  equations  of 


,aB 


(1.2.4. 5) 


0 


T 


where  the  lateral  shear  contribution  has  been  omitted.  The  semi¬ 
colon  indicates  covariant  differentiation  with  respect  to  the  Lagrangian 
co-ordinates  on  the  middle  surface  of  the  deformed  shell.  Equation 
1.2. 4. 5  is  not  so  easily  forced  to  be  an  identity  by  a  stress  function. 
Suppose  one  generalizes  Kempner’s  solution  of 


to 


where  the  comma  indicates  partial  differentiation  and  e  is  the 

a6 

permutation  symbol,  while  e  is  the  permutation  tensor.  Substitution 


into  Eq.  1.2. 4. 5  yields,  for  example, 


F  -  F 

; 2 ; 2 ; 1  ;2;1;2 


(1. 2.4.6) 


where  g  is  the  determinant  of  the  middle  surface  metric  tensor  co¬ 
efficients.  The  right-hand  side  of  Eq.  1.2. 4.6  is  not,  in  general,  zero. 
Novozhilov  (22)  gives  an  extended  discussion  detailing  when  this  equation 
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may  be  assumed  to  be  sufficiently  close  to  zero  in  small-deflection 
theory.  Langhaar  (23)  appears  to  have  found  that  Eq.  1.2.4. 5  becomes 


an  identity  for 

maS 


u<$  By  -n  ,  aB  „ 
=  £  e  *SjY  +  g  G 


where 

G  =  kF  - 

The  symbol  <  represents  the  Gaussian  curvature  of  the  deformed  middle 
surface,,  This  formulation  represents  considerable  complication  over 
Kempner 's  assumed  form.  However,  Langhaar  and  Boresi  (24)  argue  that 
Eq.  1.2. 4.5  may  readily  be  replaced  by 


since  the  part  neglected  involves  the  Christoffel  symbols.  The  Christ- 
offel  symbols,  being  functions  of  §ag»  must  approximate  their  original 
zero  values  if  the  deformation  is  largely  inextensional . 


Kempner  obtained  his  w  format  by  iterating  the  classical 

solution  in  the  following  manner.  Substitute 

(ttx  Try  ,  .  \ 

COS  Tj—  COS  Tp-  +  dl  I 

x  y 

in  Eq.  1.2.4. 4  to  obtain  Fi ,  then  use  Fj  and  wi  on  the  right-hand 
side  (i.e.,  the  nonlinear  side)  of  Eq.  1.2. 4.1  to  obtain 


W2  =  at  (  cos 

TTX 

1 

X 

TTV 

COS  Tp-  + 

2  TTX 

b  cos  i  +  c 

X 

2TTy 

cos  - — +  d 

V 

y 

y 

TTX 

+  e  cos  - —  cos 

3tt-v 

1 

-  +  g  COS 

3ttx  ry  \ 

1  COS  1 

(1.2. 4. 7) 

X 

y 

X  y 

Kempner  did  not  retain  the  last  two  terms;  he  obtained  the  final  form 
of  F  from  Eq.  1. 2.4.4  and  varied  the  potential  energy  with  respect  to 


* 


-28- 


bs  c,  ys  n  =  n2(t/R)9  and 


0) 


an 


While  studying  the  literature,  the  author  compared  Kempner*s 
format  with  the  Fourier  series  expansion  of  an  ideal  buckle  given  by 
Cox  (25)  to  be 

2  (  cos  cos  ^  +  ^(1  ~  cos  2m6  -  cos  2mc|>)j 

m  '  *  - 


W  = 


4lh 


where 


1  =  inclination  of  facet  to  axis  of  cylinder, 

h  =  axial  facet  length, 
b  =  half  width  of  facet, 
x  =  axial  co-ordinate, 


y 

e 


circumferential  co-ordinate, 
ttx  _  2ttR  yx  _  ynx 


2b  R 


R 


where  y  =  — 


and 


Then 


n  =  number  of  buckles. 


sin  y  — 


R(1  -  cos  r/n) 


ja  jt 
2  n 


so  that  the  coefficient  in  w  becomes 

4¥h  _4_  _yn_  b_  _  2R 

tt2  tt2  2n  y  n^" 

This  indicates  that  it  would  probably  be  better  to  normalize  Eq.  1.2 .4. 7 
with  respect  to  R/n2  instead  of  t.  This  normalization  was  first 
introduced  by  Yoshimura  (15)  and  more  recently  by  Hoff,  Madsen,  and 
Mayers  (12s  1966). 


Von  Kerman  and  Tsien  (8)  discussed  the  two  loading  methods 


of  prescribed  a  and  prescribed  e  (see  Fig0  12) .  Kempner  noted  that 
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u 


Kempne r '  s  stable  range  for 


prescribed  e 


Prescribed  a 


aR 

Et 


eR/ 1 

Figo  12  Loading  Methods 


the  parameters  (b,  c,  y,  r) ,  oo)  were  the  same  for  the  two  loading 
conditions  in  the  post-buckling  range,  but  that  the  stability  of  the 
state  depended  on  the  method  of  loading  (see  Fig.  12) .  He  obtained 
a  minimum  post-buckling  load  of  0.182  which  occurred  at  y  =  0.362. 

This  minimum  is  independent  of  R/t,  as  have  been  the  previous  results. 
The  value  of  (ad)t  was  found  to  be 


t  (ad) 


a  positive  quantity  which  was  approximately  t  at  the  minimum  load. 
The  maximum  radial  deflection  at  the  mimimum  load  was  8t. 

1.2.5.  Almroth  (11,  1963) 

It  can  be  shown  [see  Hoff,  Madsen,  and  Mayers  (12)]  that 
the  "Ideal''  shape  of  Fig.  1  is  described  in  part  by  the  following 
radial  deflection  from  the  middle  surface  of  a  circular  cylinder: 
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~  r  a  1TT  x  .ITTy 

£  L  A„  .  cos  — —  cos  u- 

j=0  1=0  1J  Xx  \ 


where  i  +  j  is  an  even  number „ 


Almroth  made  the  advance  of  considering  more  terms  of  the 
series  and  of  noting  their  effect „  He  made  computations  for: 

CASE  Is  Aq  o  »  A20,  An,  A40,  A22 

CASE  2:  Aoo»  A2o9  An,  A40,  A6q,  A33 

CASE  3:  AqOs  A2o  ,  An,  A40  ,  A3i  ,  A22  ,  Ai  3  ,  A60,  A33 

By  use  of  the  Newton-Raphson  method  he  was  able  to  eliminate  the  cross¬ 
plots,  etc. ,  which  had  been  used  in  previous  solutions  and  thus  considered 
a  large  number  of  terms  with  comparative  ease.  By  letting  q  be  the 
generalized  co-ordinates,  then  for  prescribed  o  it  is  required  to 
solve  the  equilibrium  equations 


3V 

3a  , 


=  0 


Almroth  included  the  variable  n  among  his  generalized  co-ordinates 
These  equations  have  the  approximate  solution 

(k+D  /  n>  r /  32v  \  (k)rVav  \(k) 

(qi)  =  lqi)  J 

where  k  and  k  +  1  denote  successive  approximations.  Almroth’ s 


basic  results  are  given  in  Fig„  13. 


fe 
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Fig.  13  Almroth's  Results 

It  is  interesting  to  note  that  Almroth  feels  that  the  variation  with 
respect  to  n  may  be  termed  as  "definitely  reasonable".  Almroth 
obtained  the  approximate  experimental  value  of  0.06  (for  cR/Et) 
and  he  stated  that  additional  terms  did  not  seem  important,  but  this 
conclusion  was  not  supported  by  the  results  of  Hoff,  Madsen,  and 
Mayers  (12).  One  notes  that  the  addition  of  four  terms,  previously 
considered  unimportant  and  also  neglected  because  of  the  intractable 
computations,  has  decreased  the  post-buckling  load  to  one-third  of 
Kempner's  value.  A  comparison  of  both  Kempner's  and  Almroth's  results 
with  Thielemann's  (21)  experimental  results  yields  excellent  correspondence 
for  axial  shortening  (see  Fig.  14). 

Again,  because  of  the  fact  that  Almroth  varied  n  and  because 
of  the  neglect  of  nonlinear  u  and  v  terms  [Tsao  (26)  has  shown  that 
retention  of  these  terms  in  the  strain  energy  introduces  an  R/ t  effect] 
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a  single  load-deflection  curve  was  obtained.  The  radial  deflection  at 
the  minimum  load  was  50t„ 


lo  2  0  6<,  Hoff,  Madsen,  and  Mayers  (12,  1966) 


It  was  the  interest  of  these  authors  to  test  the  previously 
discussed  procedures  for  a  larger  number  of  terms.  The  terms  they 
found  most  important  were  Aqqs  A1  j  ,  A02»  A2q>  A22j  A13s  A31  ,  A33, 

A40  *  A4  4»  A5  5>  A6  0»  A6 6  >  A8  0s  A100>  A120s  b>  and  n- 


format 


A  major  change  from  the  previous  papers  was  their  use  of  the 


£  A .  . 


cos 


iTTX 

1 

X 


cos 


In. 

1 

y 


w 


* 
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which  agrees  with  comments  made  previously  concerning  Kempner’s 
paper.  The  significance  of  this  form  is  that  the  exact  diamond  pattern 
has  A. .  as  a  function  of  i  and  j  alone,  and  not  of  1  and  n;  so 
that  if  the  derived  shape  approximates  a  diamond  pattern,  then  the  A 
values  should  approach  the  theoretical  values  independent  of  n  and  1  « 
Hoff,  et.  al.,  found  that  the  A  values  appeared  to  converge  to  the 
ideal  values.  It  is  to  be  noted  that  Yoshimura  (15)  first  suggested 
this  form  of  normalization,  and  he  also  commented  upon  the  constancy 
of  the  coefficient  values  so  obtained. 

The  maximum  lateral  deflection  was  found  to  be  236t.  Figures 

15a,  15b,  and  15c  contain  the  major  results.  The  numbers  in  brackets 

indicate  the  number  of  A. .  coefficients  considered. 

ij 


1  2  3  4  5 

eR/t 

Fig.  15a  aR/Et  versus  eR/t 


0  1  2  3  4  5  6 

eR/ 1 

Fig.  15b  y  versus  eR/t 


0 
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n 
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0  1  2  3  4  5  6 

eR/t 

Fig.  15c  n  versus  eR/t 

It  is  clear  that  as  the  number  of  terms  retained  increases, 
oR/Et,  n ,  and  y  all  shift  downward,  as  indicated  by  the  arrows.  Hoff, 
et .  al.,  concluded  that  the  use  of  an  infinite  number  of  terms  of  the 
series  would  yield  the  result  of  an  exact  diamond  shape  for  which ; 

i.  n  =  0  and  since  n  >_  2  this  implies  t/R  =  0 

ii  o  y  =  0  which  implies  n  =  00 

and 

iii.  oR/Et  =  0  which  implies  a  =  0  ; 

in  short,  a  trivial  result. 

Since  there  is  no  logical  objection  to  using  all  of  the 
terms  (there  is,  as  noted,  an  objection  to  minimizing  with  respect 
to  n  and  y) ,  and  since  the  experiments  of  Thielemann  (21)  implied 
that  there  exists  a  non-trivial  minimum  load  in  the  large-deflection 
state,  then  the  conclusion  would  seem  to  be  that  the  method  fails. 


Further  discussion  on  this  point  shall  be  reserved  until  Chapter  4 
when  the  results  of  work  done  by  Jones  (16)  and  Hoff  and  Madsen  (27) 
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will  be  considered  and  compared  with  the  results  of  the  present 
work,, 


Hoff,  et.  al.,  commented  upon  the  relative  inaccuracy  of  the 
present  strain  tensor  for  large  deflections.  Let  us  consider  the 
example  given  by  them  for  n  equals  six.  In  such  a  case,  Donnell’s 
strain  relations,  applied  to  the  ideal  diamond  shape  (assuming  steel), 
gives  a  stress  of  300,000  p.s.i.  Yoshimura  (15)  gave  a  plot  of  the 
middle  plane  stress  distribution  for  y  =  0.7  and  n  =  9.4, 
which  revealed  axial  tension  stresses  which  were  two  thirds  the 
magnitude  of  the  axial  compressive  stresses.  Legget  and  Jones  (20,  1942) 
and  Madsen  and  Hoff  (27,  1965)  have  reported  similar  stress  distri¬ 
butions.  One  does  not  expect  the  axial  middle  plane  stresses  to  become 
tensile.  It  is  felt  that  these  results  show  that  the  von  Kdrman-Tsien 
model  is  not  applicable  to  the  post-buckling  problem.  That  Hoff,  et .  al., 
probably  did  not  have  axial  tension  stresses  must  be  attributed  to 
their  value  of  y,  which  gives  an  axial  wave  length  approximately  six 
times  that  indicated  by  experiment. 

1.2.7.  Uemura  (13,  1963) 

It  was  previously  noted  that  buckles  usually  occur  over 
only  a  portion  of  the  cylinder.  Uemura  considered  a  large-deflection 
approximation  to  the  "Tier"  configuration  of  Fig.  2a,  and  he  made  use 
of  an  asymptotic  unperiodic  w  displacement  function.  His  procedure, 
which  is  involved  algebraically  and  not  sufficiently  explained  with 
respect  to  method  and  approximations,  will  not  be  discussed,  save  to 


-36- 


mention  that  he  eventually  minimized  with  respect  to  four  variables, 

one  of  which  was  n.  Uemura  found  that  nearly  square,  damped  buckles 

could  occur  at  a  post-buckling  load  of  approximately  0.5a  .  The 

cl 

results  were  dependent  upon  /rF/L  and  so  exhibited  an  L/R  effect. 
The  solution  did  not  maintain  complete  geometric  continuity.  Yoshimura 
(15),  who  presented  a  qualitative  analysis  of  local  buckling,  was  also 
unable  to  maintain  geometric  continuity  in  his  solution. 

1 o  2 . 8 .  Pogorelov  (14,  1962) 

In  a  series  of  papers,  Pogorelov  considered  thin-shell 
buckling  based  upon  a  neglecting  of  the  membrane  energy.  The  buckled 
form  considered  was  approximately  isometric  to  a  cylinder  and  was 
constructed  by  way  of  a  mapping  procedure  through  a  2n  regular 
prism.  The  term  "approximate"  is  used  because,  like  Yoshimura  (15) 
and  Uemura  (13)  ,  Pogorelov  was  unable  to  maintain  complete  geometric 
continuity.  Pogorelov' s  arguments  are  largely  geometric  in  nature  and 
are  strictly  qualitative. 

Pogorelov  assumed  that  the  post-buckling  wave  was  square 

(i.e.,  y  =  1)  and  that  n  was  determined  from  the  solution  for  the 

minimum  classical  neutral  load  under  the  boundary  conditions: 

i.  w  =  0  , 

ii.  M  =  0  , 

xx 

iii.  v  =  0  , 

and 


iv. 


buckling  occurs  under  constant  load. 
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Under  the  assumption  of  a  square  wave,  the  classical  relation  between 
y  and  n  gives 

1 

n(t/R)T  ^  0.9 

Thus  n  and  y  were  assumed  by  Pogorelov  to  be  established  before  the 
post-buckling  era.  By  way  of  contrast,  this  study  will  use  n  and  y 
as  the  fundamental  unknowns  and  the  results  indicate  that,  if  one 
wishes  to  determine  a  stable  post-buckling  configuration,  there  is 
good  reason  to  assume  n  to  be  initially  established. 

Pogorelov* s  analysis  was  based  on  the  assumption  that  the 
undeformed  cylinder  was  the  virgin  state  of  the  material,  and  covered 
three  distinct  stages: 

1.  Classical  analysis:  The  classical  analysis  determined 
n  under  the  assumption  of  a  square  buckle  pattern. 

2.  Unstable  post-buckling  era:  Using  y  =  1,  Pogorelov 

found 

aR/Et  =  0.138  , 

but  the  solution  was  unstable. 

3.  Stable  post-buckling  era:  Using  y>0.91  (the  exact 
value  of  y  was  not  determined) ,  Pogorelov  obtained 

aR/Et  =  0.15 

Since  the  classical  n  exceeds  the  experimental,  a  better  value  of 
n  might  possibly  be  found  from  alternate  solutions  to  the  neutral 
equilibrium  equations  [see  Ohira  (28),  Hoff  (29),  Hoff  and  Rehfield  (30)]. 
However,  the  highspeed  photographic  experiments  of  Almroth,  Holmes, 
and  Brush  (31)  have  shown  that  the  number  of  initial  buckles  decreases 
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by  half  when  the  stable  state  is  reached,  therefore  indicating  that  the 
classical  n  does  not  apply  to  the  post-buckling  problem. 

It  is  clear  that  Pogorelov’s  procedure  differs  radically 
from  that  of  the  previously  discussed  works;  indeed,  this  author  has 
not  read  a  major  American  work  which  even  referenced  his  contributions. 
Pogorelov’s  work  is  of  interest  to  this  author  because  Pogorelov's 
methods  are  comparable  to  his  own  and  are,  from  the  outset,  designed 
as  qualitative  approaches. 

1.2.9.  Closure  to  the  Literature  Discussion 

Thielemann  (32,  1960),  in  a  survey  of  thin-cylinder  buckling, 
suggested  that  further  advances  could  possibly  consist  of  the  following 
a  more  accurate  analysis  of  initial  imperfections  with  a 
view  to  reducing  the  upper  critical  load; 

the  use  of  a  w  function  involving  more  terms  and  a  deter¬ 
mination  of  the  most  important  terms; 

and  the  use  of  asymptotic  unperiodic  displacement  functions 
to  cover  local  buckling. 

Hoff,  et.  al. ,  (12,  1966)  considered  the  use  of  a  more 
complete  w  function  in  detail  and  Uemura  has  initiated  a  study  of 
the  use  of  asymptotic  unperiodic  displacement  functions.  The  results 
of  Hoff,  et.  al.,  are  most  surprising  and  indicate  that  the  current 
procedures  lead  to  a  trivial  result.  The  author  is  not  aware  of  an 
initial  imperfection  analysis,  using  large-deflection  theory,  which 
differs  significantly  from  that  of  Donnell  and  Wan  (9).  The  upper 
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critical  load  has  received  much  consideration  by  means  of  the  neutral 
equilibrium  approach.  If  these  equations  are  solved  subject  to  the 
boundary  conditions 

u  =  0,  v  =  0,  M  =  0, 

xx 

and  if  the  buckling  is  assumed  to  occur  under  constant  load,  then  a  , 

cl 

is  obtained.  Hoff  (29),  Hoff  and  Rehfield  (30),  and  Ohira  (28),  have 
shown  that  enforcing  other  reasonable  boundary  conditions  leads  to 
approximately  a  fifty  percent  reduction  in  the  neutral  load.  If  one 
notes  that  the  neutral  partial  differential  equations  are  linear,  and 
that  the  post-buckling  total  deflection  equations  are  highly  nonlinear, 
then  it  would  seem  that  the  more  accurate  solution  of  the  linear 
system  is  the  easier  task.  From  a  structural  point  of  view  the  non¬ 
linear  analysis  would  seem  to  be  the  most  desirable.  Mathematically 
and  physically  both  problems  are  of  interest. 

In  conclusion,  the  preceding  discussion  of  the  literature 
has  made  evident  the  following: 

1.  A  comparison  of  the  theoretical  results  with  experimental 
evidence  has  shown  that  the  theoretical  solutions  must  be  judged  as 
qualitative  and  not  quantitative  as  some  authors  imply. 

2.  Despite  the  passage  of  nearly  thirty  years  since  the 
problem  was  first  defined  by  von  Kdrman  and  Tsien,  there  is  still  no 
generally  accepted  method  by  which  to  define  the  values  of  p  and  n. 

3.  The  ability  of  the  von  Karman-Tsien  model  to  represent 
the  large-deflection  post-buckling  state  is  to  be  questioned.  Von 
Karmdn  and  Tsien  (8,  1941)  and  Donnell  and  Wan  (9,  1950)  have  pointed 
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this  out,  yet  many  authors  continue  to  apply  the  model  without,  it 
seems,  sufficient  justification. 

4.  The  works  of  Legget  and  Jones  (20,  1942),  Yoshimura 
(15,  1951),  and  of  Madsen  and  Hoff  (27,  1965)  have  shown  that  the 
von  Karman-Tsien  model  may  yield  tensile  values  of  T  .  One  does 

XX 

not  expect  tensile  values  of  T  .  This  indicates  the  need  to  include 

xx 

the  T  distribution  in  published  reports  (as  well  as  the  value  of 
xx 

aR/Et)  in  order  that  the  solution  may  be  properly  assessed. 

5.  The  use  of  modern  computation  techniques  has  allowed 
the  consideration  of  many  variables;  however,  the  conclusions  reached 
cannot  be  regarded  as  final. 

1.2.10.  Purpose  of  the  Thesis 

It  is  the  purpose  of  this  study  to  further  consider  the 
post-buckling  behavior  of  an  axially  compressed  thin  circular  cylind¬ 
rical  shell.  For  reasons  to  be  explained,  the  analysis  will  largely 
neglect  the  middle  plane  deformation  and  will  consist  of  "patching" 
the  deformed  surface  together  as  suggested  by  Kennedy  (33).  In  his 
treatment  of  the  middle  plane  deformation,  the  author  of  this  thesis 
recognized  the  difficulty  of  obtaining  a  large  deflection  mapping  which 
would  map  the  middle  surface  almost  inextensionally ;  thus  it  was  the 
aim  of  his  procedure  to  take  note  of  whether  or  not  the  y  value 
so  obtained  would  be  closer  to  the  experimental  results  than  is  the 

U 


value  of  Hoff,  et,  al.  ,  (12). 


, 
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The  w  deflection  function  of  Hoff,  et.  al.,  (12)  will  not 
be  assumed;  rather,  the  large  deflection  form  will  be  approximated  by 
zones  of  bending  (described  in  closed  form  as  spiral  curves) ,  plus  a 
flat  portion.  The  undeformed  cylinder  will  be  assumed  to  be  developed 
from  a  flat  plate.  It  seemed  desirable  to  express  the  bending  zone  in 
closed  form  with  the  hope  that  the  number  of  parameters  could  be  reduced 
and  the  analysis  simplified.  Unfortunately,  it  was  not  possible  to 
maintain  complete  geometric  continuity  and  this  probably  resulted  in 
a  less  stiff  structure.  As  mentioned  in  the  Introduction,  the  model 
presented  employed  three  parameters  and  yielded  results  comparable 
to  present  theory. 


CHAPTER  2 


GEOMETRY  OF  THE  LARGE  DEFLECTION  MODEL 


It  is  the  purpose  of  this  chapter  to  establish  an  approxi¬ 
mate  large  deflection  geometry  for  the  nost-buckling  configuration  of 
a  thin  circular  cylindrical  shell  under  axial  load. 


The  entire  cylinder  is  assumed  to  have  buckled  into  a  periodic 
diamond  pattern  and  no  end  effect  is  considered.  The  unbuckled  shell 
will  be  assumed  to  be  developed  from  a  flat  plate.  Bending  will  be 
considered  to  be  the  primary  mode  of  deformation;  the  zones  of  bending 
are  indicated  below  in  the  diagram  of  the  ideal  buckle  pattern  (see 
Fig.  16).  A  typical  cross  section  is  shown  in  Fig.  17. 


Rib 


Flat 


This  cross 
section  is 
shown  in 
Fig.  1 


Flat 


Zones  of 
Bending 


Fig.  16  Idealized  Shape 
-42- 


-43- 


Flg.  17  Cross  Section  of  Fig.  16 


It  is  clear  that  the  ideal  portion  of  Fig.  17  is  valid  only 
if  t  3  0  (i.e.,  no  bending  stiffness);  in  the  actual  approximation 

a  transition  curve  of  finite  curvature  must  be  provided  between  adjacent 
flats.  In  the  following  analysis,  the  transition  curve  will  be  assumed 
to  be  a  spiral. 

2.1.  Discussion  of  the  Spiral  Curve 

The  spiral  curve  possesses  the  property  that 
rl  *  RL  .  (2.1.1) 

With  reference  to  Fig.  18,  let 

*  f(4>) 


1 
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The  differential  equation  of  the  curve  is 

dl  RL 
r  d<j>  1 

and  integration  yields 

l2 

—  =  RL<|)  +  constant  , 

which,  if  1  =  0  when  <f>  =  0,  becomes 

l2  =  2RLcf>  .  (2.1.2) 

From  Eqs.  2.1.2  and  2.1.1  it  is  evident  that  at  <p  =  0,  r  =  00 

and  that  at  4)  =  A,  r  =  R;  hence,  the  spiral  is  useful  as  a  trans¬ 
ition  curve  between  a  straight  line  (the  flat)  and  a  curve  of  prescribed 
radius  (at  the  rib) . 


The  cartesian  co-ordinates,  expressed  parametrically  in 


terms  of  <j>,  are  obtained  as  follows.  From  Eq.  2.1.2 

1 


dl  = 


A<p 


d<|> 


(2.1.3) 


where 


Since 


then 


A  =  (2RL)2  =  2RA2 


dx 

dl 


=  COS<J) 


and 


dl 


=  sin<J) 


dx  = 


M 
2 


,1_ 

2 


cosc|)d<j) 


and  dy  = 


A  b 


.1_ 

J 


(2.1.4) 


sin<J>d<j> 


When  the  sine  and  cosine  functions  are  replaced  by  their  Taylor's 
series  expansion  about  <j>  =  0  and  the  resulting  expressions  are 

integrated,  the  x  and  y  co-ordinates  are  found  to  be 


x  =  A<f> 


!  _  JL.  H  _  .  j>.6-  +  ) 

1  5.2!  9.4!  13.6!  *  *  */ 


(2.1.5a) 
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and 


y  = 


1 


i  _ 

3  7.3! 


+ 


4. 


11.5! 


i! _ 

15.7! 


C 


The  distance  to  any  point  on  the  curve  from  the  origin 


given  by 


1 


» 


c 


and  the  total  length  of  the  spiral  is 
L  =  2RA 

The  angle  A,  as  shown  in  Fig,  18,  represents  the  total  angle  of 
turning  of  the  curve. 


2.2.  General  Procedure  to  be  Followed  in  Developing  the  Mapping 


.1.5b) 

is 

.1.6) 


Functions 


Figure  19  explains  some  necessary  geometry.  The  diamond 
pattern  is  to  be  mapped  onto  the  space  frame  as  shown  in  Fig.  19, 
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where  R  equals  the  initial  cylinder  radius  and  MR  equals  the  rib 
radius . 


Z 


Fig.  19  Space  Frame 


A  typical  diamond  of  the  buckled  shell  will  appear  as  shown 
in  Fig.  20.  Figure  20  serves  to  define  some  of  the  necessary  notation. 


Fig.  20  Diagram  of  State  1 
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If  n  is  the  number  of  buckles,  then 

2ttR 

n 


B0  = 


and 


Lq  =  BBq 


where  the  relation  between  y  and  B  is 


B  = 

2y 


(2.2.1) 


(2.2.2) 


As  indicated  in  Fig.  20,  the  zones  of  bending  in  the  edge 
and  middle  spirals  have  been  arbitrarily  taken  of  equal  width  d; 
this  simplified  the  computations  and  forced  the  geometric  constraint 
discussed  below. 


A  series  of  mappings  using  Lagrangian  co-ordinates  (x,y), 
initially  the  orthogonal  cartesian  co-ordinates  of  the  middle  plane 
of  the  flat  plate,  state  1,  will  now  be  developed  which  will  map  the 
flat  plate  (shown  in  Fig.  20)  onto  the  space  frame  (shown  in  Fig.  19). 
The  middle  surface  will  be  mapped  inextensionally  in  the  x  direction; 
the  Lagrangian  x  co-ordinates  of  the  edge  spiral  at  the  rib  and  at 
the  transition  are  x^  and  xT.  The  mapping  functions  will  not  be 
periodic;  to  this  end,  the  geometric  constraints  of  c  0  and 

c  4=  00  must  be  imposed.  (See  Eq.  2.5.1  for  the  definition  of  c.) 

The  geometric  constraint  on  c  accounts  for  the  requirement  that  the 
plate  was  to  be  mapped  into  a  closed  cylinder  which  should  buckle 
into  a  periodic  closed  shape  with  continuous  tangents.  The  net 
result  should  then  approximate  a  diamond-shaped  buckle.  This  mapping 
is  normally  effected  by  the  use  of  a  displacement  U  from  the  initial 
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cylinder;  however,  this  procedure  was  abandoned  since  the  easiest  way 
to  get  U  seemed  to  be  the  method  used.  One  could,  for  instance, 
specify  a  kinematically  admissable  w  displacement  function  with 
respect  to  a  local  reference  frame  based  on  the  middle  surface  of  the 
undeformed  cylinder,  and  then  obtain  v  by  the  constraint  that 

eH  =  ° 

where  e^  is  defined  by  Eq.  2.9.6.  Assuming  that  u  has  a  small 
effect  in  this  equation,  then  the  equation  may  be  considered  as  a 
nonlinear  partial  differential  equation  defining  v.  An  unsuccessful 
attempt  was  made  to  use  this  method.  The  differential  equation  was 
integrated  by  a  method  of  successive  approximations;  the  procedure 
was  involved,  difficult,  and  made  the  consideration  of  many  terms 
necessary.  However,  it  was  found  that  the  procedure  used  in  section  2.3, 
etc.,  converted  this  complicated  problem  to  a  problem  which  was  readily 
and  easily  solved.  The  mapping  technique  of  section  2.3  obviated  the 
need  to  solve  a  nonlinear  differential  equation  in  order  to  enforce 
the  inextensional  constraint  in  the  x  direction. 

The  basic  co-ordinate  frame  is  orthogonal  cartesian  (X,  Y,  Z) 
and  is  oriented  as  shown  in  Fig.  19.  The  various  mappings  will  be 
called  "states"  and  the  corresponding  basic  co-ordinates  will  be  suit¬ 
ably  subscripted  [say,  (X  ,  Y1 ,  Z^) ] .  The  additional  subscripts  of  R 
and  T  will  denote  rib  and  transition  values  respectively;  for 

example,  (X^^,  Y.^,  ^1R^ ’  ^IT*  ^1T*  ^1T^  * 
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2.3.  Mapping  of  the  Edge  Spiral 


and 


and 


The 


transition  equations  for  state  1  are: 


=  y 


The  transformation  equations  for  state  2  are: 

^2  ~  x  * 

Y2  =  Ny  , 


(2.3.1a) 

(2.3.1b) 


(2.3.1c) 


(2.3.2a) 

(2.3.2b) 


(2.3.2c) 


State  3  is  shown  in  Fig.  21.  Assign  to  the  co-ordinates 

(X^,  Z3R) •  equations  of  the  spiral  in  terms  of  (X1,  Z1), 

see  Fig.  21,  are  known  by  Eqs.  2.1.5. 


Fig.  21  State  3  at  y  =  0 
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If  the  middle  line  is  mapped  inextensionally ,  then 
1  =  |x  -  x  | 

=  xT  -  X  , 


so  that 


*  =  (xT  ‘  X)2 

a"2 


(2.3.3) 


Letting 


and 


then 


X”  =  -X1 


z"  =  -z1 


where 


X, 


x3R  +  (X"  -  X-  )  , 


(2.3.4a) 


and 


Y  =  Y 
3  2 


Z3  ■  Z3R+(Z"-ZV’ 


(2.3.4c) 


where 


X"  =  -AA2  (1  -  T-vr  + 


A 


R 


5.2!  9.4! 


13 


.6!  +  “•) 


(2.3.5a) 


and 


Z" 


R 


1 

-AA2 


/A  _  A3  A^ _ A^ _  A 

\3  "  7.3!  11.5!  15.7!  ***  J 


(2.3.5b) 


A  =  - 

n 


(2.3.6) 


Further , 


X 


ii  _ 


-(x  -x)  (  !  .  (XT  -x)l*  ,  (xT  -  x)6  (xT  -  X)12  + 

T  '  5 . 2 ! A4  9 . 4 ! Aa  13 . 6! A1 z 


(2.3.7a) 
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and 


Z"  = 


-(xT  -  x)  |- 


(XT  -  x)2  (XT  -  x)6,  (XT  -  x)1Q 


3AZ 


7 . 3 !  Ab  +  11~5 !  A1  b  tTTTTnr  +• 


The  co-ordinates  of  xg  are  given  by 

X  =  -  MR  ^  sin  A  , 

O  iv  L  q 


and 


Y3R  '  NY  • 


Z_  =  MR  -  MR  ^ 


3R 


L0 


(1  -  cos  A) 


(XT  -  x)14 
15 .  7 !  AJ 


(2.3.7b) 

(2.3.8a) 

(2.3.8b) 

(2.3.8c) 


2.4.  Mapping  of  the  Flat 

The  co-ordinates  of  the  transition  are  obtained  by  evaluating 
Eqs.  2.3.4  at  X"  =  Z"  =  0  .  Hence 


X3T  XT  ’ 


and 


3T 


'3T 


=  Ny 


=  Z 


3R 


ZR 


The  co-ordinates  of  a  general  point  on  the  flat  become, 
X3  -  x  . 


Y3  =  Ny 


and 


Z3  Z3T 


(2.4.1a) 


(2.4.1b) 


(2.4.1c) 


(2.4.2a) 


(2.4.2b) 


(2.4. 2c) 


2.5.  Determination  of  the  Independent  Parameters 


The  radius  of  the  edge  spiral  at  the  rib  is  defined  by 


R  =  cR 
e 


(2.5.1) 
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The  parameters  introduced  so  far  are  n,  N,  M,  B,  and  c. 
However,  they  are  not  all  independent,  as  will  now  be  shown.  The 
mapping  of  the  edge  spiral  and  flat  has  been  under  the  constraint 
that  the  middle  line  was  neutral;  hence,  at  y  =  0,  there  must  hold 

the  relation 


T  ttR 

L  -  XT  "  7T 


(2.5.2) 


From  continuity  in  the  X  direction  it  is  also  required  that 


X3R  XT  +  X"r 


(2.5.3) 


Substituting  for  x^,  from  Eq.  2.5.3  into  Eq  0  2.5.2  results  in 


L  -  (X3R  -  XV  '  f 


(2.5.4) 


Replacing  X  by  Eq.  2.3.8a  and  X"  by  Eq.  2.3.5a  and,  noting  that 
JK  R 


L  =  2cRA,  reduces  Eq.  2.5.4  to 


^  =  c2RA 

n 


-  [- 


MR  sinA  +  AA2  (  1  - 


A" 


A 


5.2!  9.4 !  13 


—  +  11 
.61  +>--j  J 


which,  upon  replacing  A  by 

1_ 

A  =  2cRA2 


(2.5.5) 


becomes 


c  = 


A  - 


2A  -  2A  1  - 


M  sinA 

2 - 


A 


5.2, 


+ 


A 


A 


9.4!  13.6! 


+ 


2 


(2.5.6) 


Therefore,  c,  M,  and  n  are  not  independent.  As  indicated  in  section  2.2 
the  parameter  c  is  subject  to  a  geometric  constraint;  hence,  Eq .  2.5.6 
serves  to  define  M  =  M(c,n). 


2.6.  Calculation  of  the  Base  Vectors  for  the  Edge  Spiral  and  for  the  Flat 

The  position  vector  of  a  typical  point  on  the  middle  surface 
of  the  spiral  is 

R  =  X3  1  +  Y3±  +  Z3K  ,  (2.6.1) 
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where  the  components  are  given  by  Eqs.  2.3.4. 


which , 


The  equation  of  the  x^  line  is 

B0  -  2L  ^ 

~  )  L° 

solving  for  x^,  becomes 

xT  =  [y  -  (1  -  2c)  2ABR] 


y  = 


2L( 

b7 


XT  + 


Defining 


y  =  (1  -  2c)  2ABR 


then  is  given  by 


xn 


=  (y_  y) 


2B 


This  is  the  value  of  x^  to  be  substituted  into  Eq.  2.6.1, 


(2.6.2) 


(2.6.3) 

(2.6.4) 


By  definition,  the  base  vector  on  the  middle  surface  and  in 


the  x  direction  is 

7  -  il 

Go1  "  9x 


(2.6.5) 


Performing  the  computations  indicated  by  Eq .  2.6.5  results  in 

(x  -  x)4  (x  -  x)8  (x  -  x)12 

— -  +  - -  -  -  +.  .  . 

2IA4  4! A8  6|A12 

/  (x  -  x) 2  (x  -  x)6  (x  -  x)lu  (x  -  x) 1 4  V 

+  — -  -  -  +  — -  -  —  +.  ..  ]  K 

\  p2.  3  «  a6  5  »  a1  0  ?  »  a1  i  / 

By  use  of  Eq.  2.3.3  (i.e.,  <p  =  (x^  -  x)2/A2),  this  equation  simplifies 

to 

G01  =  cos  <J>  i.  +  0j_  +  sin  tp  K  .  (2.6.6a) 

Thus,  while  the  (X,  Y,  Z)  co-ordinates  of  the  spiral  are  not  expressible 
in  closed  form,  the  constraint  that  the  middle  line  is  mapped  inexten- 
sionally  allows  the  base  vector  to  go  into  closed  form.  This  is  a 


1  +  °1 
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noteworthy  simplification,  for  otherwise,  the  computation  of  the 
metric  tensor  components  would  involve  the  multiplication  of  infinite 
sums . 


Proceeding  in  a  similar  manner,  one  may  derive  that 


-  _  3R 

G°2  “  37 


MR  .  .  cosd) 
—  slnA 


\  MR  ,,  A  \  sin<}>  \  Tr 

)  i  +  NJ.+  (1  -  cosA)  “  2B~^  )  — 

(2.6.6b) 

Analagous  computations  for  the  flat  yield 

Gqi  =  U  +  pi  +  OK  ,  (2.6.7a) 


and 

_  MV 

G02  =  Oi  +  Nj_  4-  —  (1  -  cosA)  K  .  (2.6.7b) 

L0 

The  continuity  of  the  base  vectors  may  be  checked  by  comparing  Eqs.  2.6.6, 
evaluated  at  <f>  =  0,  with  Eqs.  2.6.7.  The  former  yield 

Gqi  =  li  +  0 j_  +  OK  , 


and 


'02 


/  MR  .  .  1  \  .  MR  A  \  TT 

- —  sm  A  -  —  I  l  +  N i  +  - —  (1-  cosA)  K 

\  Eq  ■‘-E  /  Eq 


Continuity  clearly  requires  that 

MR  .  .  i  n 

17  Sln  4  -  2B  "  ° 


or 

M  sin  A  -  A  =  0  .  (2.6.8) 

If  Eq.  2.6.8  is  obeyed,  then  Eq .  2.5.6  yields  c  =  0,  i.e.,  the 

perfect  diamond  shape.  Therefore,  while  the  surface  has  been  mapped  in 
a  continuous  fashion,  the  material  in  the  spiral  has  been  sheared.  The 
shearing  is  exactly  similar  to  that  which  would  be  required  to  cause 
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a  piece  of  wax  paper  (Fig.  22)  to  roll  onto  the  cylinder  without 
advancing  down  the  axis  of  the  cylinder.  This  deficiency  was  not 
easily  corrected,  without  introducing  other  discontinuities,  and 
resulted  in  the  need  to  neglect  the  middle  plane  shear  energy. 

Equation  2.6,8  indicates  that  a  better  definition  of  M  would  be 

I 

M  =  kM0  ,  (2.6.9) 


where 


M0  = 


sinA 


The  magnitude  of  k  is  a  measure  of  how  close  the  exact  diamond 
configuration  is  approached;  k  =  1  is  required  for  the  perfect 

diamond  and  k  <  1  is  required  for  an  approximation  to  the  exact 
diamond  pattern. 


Fig.  22  Explanation  of  Shear  Discontinuity 


2.7.  Transformation  Equations  for  the  Middle  Spiral 

Figure  23  shows  the  geometry  of  the  middle  spiral.  It  was 
previously  decided  to  let  the  middle  spiral  length  be  d  (see  Fig.  20); 
hence , 

-  2V  • 


d 


(2.7.1) 
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Fig.  23  Geometry  for  Middle  Spiral 


is  the  radius  of  the  middle  spiral.  States  1  and  2  are  as  previously 
given  and  the  final  mapping  is  obtained  as  follows.  Introduce  local 
intrinsic  co-ordinates  (x,  z) ,  where 

z  =  -(d  -  y)  (2.7.2) 

State  2  will  be  mapped  maintaining  its  length.  This  is  in  error,  since 
the  curvature  introduced  must  change  the  length,  but  must  be  accepted 
as  an  approximation  which  will  be  valid  for  small  f  .  Thus,  the  co¬ 
ordinate  z  corresponds  to  a  length  of 

1  =  -Nz  (2.7.3) 


which,  using  Eq.  2.1.6,  is  also  expressible  as 
1  =  a<f>2  , 


where 

a  =  2V' 

The  angle  <f>  is  therefore  given  by 


1_ 

,2 


<f> 


N2  (— z)2 
a2 


(2.7.4) 


(2.7.5) 


(2.7.6) 
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From  Eq „  2.1.5, 

Y  =  -a<f> 


1  - 


4. 


5.2!  9.4!  13.6! 


r  +. 


and 


1 


Z  =  at}) 


7 


2  I  A  _ 

3  7.3!  11.5!  15.7! 


T  +  • 


The  transformation  equations  from  (Y,  Z)  to  (Y3  ,  Z3  )  are 


and 


Thus 


Y3  =  Y  cos  -  Z  sin  1 


Z3  =  Y  sin  Z  cos  f 


X,  =  x 


Y3  =  Y3  +  A (y) 


and 


Z  3  =  z3t  +  A(z)  +  Z  3 

The  position  vector  of  a  point  (x,  z)  on  the  spiral  is 


R  =  xi  +  (Y§  4-  A  (y) )  ±  +  (Z3t  +  A(z)  +  Z§)  K 


from  which 


9R 


Go  1  =  ^  ~  +  °i  + 


and 


9R 


9R  9z 


^°2  9y  9z  9y 

=  Oi^  +  N  (cos  1  cos  <}>  +  sin  V  sin  t}>)  j_ 

+  N  (sin  ¥  cos  cf)  -  cos  1  sin  <p)  K 

2.8.  Position  Vector  and  Base  Vectors  of  a  General  Point 


Let  z  be  a  co-ordinate  perpendicular  to  the  middle 


(2.7.7a) 

(2.7.7b) 

(2.7.8a) 

(2.7.8b) 

(2.7.9a) 

(2.7.9b) 

(2.7.9c) 

(2.7.10) 

(2.7.11a) 


(2.7.11b) 

surface 


of  the  initial  flat  plate.  Its  range  is  -t/2  <_  z  £  +  t/2,  where  t  is 
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the  thickness.  As  is  usual  in  shell  theory,  z  will  be  assumed  normal 
to  the  deformed  surface  [hence,  (x,  y,  z)  are  Gaussian  co-ordinates 
based  on  the  middle  surface]  and  unstrained. 


The  position  vector  of  a  typical  noint  (x,  y,  z)  is 
R  =  R0  (x,  y)  +  z  G03  (x,  y)  , 


(2.8.1) 


where  Rq  is  the  position  vector  of  (x,  y,  0).  Gq3  is  the  unit  normal 

to  the  deformed  middle  surface  and  is  defined  by 

g03  ^  Gqi  X  Go 2  •  (2.8.2) 

Evaluating  Eq.  2.8.2  for  the  edge  spiral  one  obtains, 

—  MR  .  .  •  A  MR  As  ,  )  . 

Gq3  =  -N  sm  cf>  _i  +  ^  —  sm  A  sm  cj>  -  —  (1  -  cos  A)  cos  cf>  J 


+  N  cos  4>  K 

Similar  computations  for  the  flat  yield 


(2.8.3) 


'0  3  =  °i  +  (  -  (1  -  cos  A)]  1  +  NR 


')  * 


(2.8.4) 


and  for  the  middle  spiral 

Go  3  =  Oi  +  N  (cos  V  sin  <p  -  sin  T  cos  <j>)  j_ 

4-  N  (cos  T  cos  +  sin  T  sin  <j>)  K  .  (2,8.5) 

Using  Eq.  2.8.1  the  base  vectors  at  a  general  point,  (x,  y,  z) ,  are 
calculated  to  be 

3D  _  _ 

,  (2.8.6a) 


—  9R  —  9Gp  3 

Gi  =  t—  =  Gqi  +  z  — 

1  9x  9x 


-  9R 

°2  =  a7  =  go2  +  2 


9G 


H3. 


9y 


(2.8.6b) 


and 


9R 

9z 


03 


(2.8.6c) 


• 
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The  base  vectors  for  the  edge  spiral  become: 

1 


cos  <J>  +  z 


2N  r 


COS  <f) 


1 


1 

2MR  .  .  T  2 MR  aN  AJ  .  \ 

s m  A  <p^  cos  <j)  -  - — —  (1  -  cos  A)  (p^  sm  <f> 


L0A 

L,-„  *  4.  ,  (  2££ 


+  sin  <j)  +  z 


A 


sin  (p 


L0A 

K 


1 

(2.8.7a) 


1 


G2  = 


MR  .  A  coscj) 

-  S.H  A  -  — 


+  z 


+ 


N  +  z 
MR 


MR 


BL0A 


sin  A  <f> z  cos  <p  + 


N4>^ 

BA 

MR 


COS  <j> 


1_ 

T 


(1  -  cosA)  cpz  sin  <() 
BL0A  i 


L0  (1  '  COS  A)  -  ^  ) 


+  z 


M 

BA 


sin  tp 


K 


1 

,  (2 


and 


G3  =  -  N  sin  <J)  _i  + 


MR  ...  .  MR  , ,  .  ,  ,  1  . 

- —  sm  A  sin  <p  -  - —  (1  -  cos  A)  cos  <J>  j 

Lo  Lo 


+  N  cos  <f>  K 

The  base  vectors  for  the  flat  become: 

G^  =  li  +  0j_  +  OK  , 

MR 


and 


G2  =  Oi  +  Ni  +  —  (1  -  cos  A)  K 

L0 


MR 


(2.8.7c) 

(2.8.8a) 

(2.8.8b) 


G3  =  Oi  -  - —  (1  -  cos  A)  i  +  NK 
“  Lo 

The  base  vectors  for  the  middle  spiral  become: 


(2.8.8c) 


Gi 

G2 


and 


1 1  +  0 ±  +  OK 

Oi  +  (cos  ¥  cos  <f>  +  sin  ¥  sin  )  Dj_ 

+  (-cos  ¥  sin  (p  +  sin  ¥  cost{>)  DK 

Oi^  +  N  (cos  ¥  sin  <f>  -  sin  ¥  cos  (f>)  J_ 
+  N  (cos  ¥  cos  <J)  +  sin  ¥  sin  <J>)  K 


(2.8.9a) 


(2.8.9b) 


.8.7b) 


(2.8.9c) 
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where 


D  = 


Nz 

N  +  z  - 

hi 


(U 


1 

T 


(2.8.10) 


2.9.  Calculation  of  Strain  Tensor  Components 


The  metric  tensor  components  are  given  by 


G.  .  =  G„  ‘  G. 

iJ  i  J 


(2.9.1) 


By  construction 


g13  ~  g23  =  0 


(2.9.2a) 


and 


G3  3  =  1 


(2.9.2b) 


Terms  containing  z2  will  be  neglected  in  G_^  .  The  remaining 


metric  tensor  components  for  the  edge  spiral  are: 
=>.  ■  ■ 


(2.9.3a) 


G22  “  2 


M 


2AB 


)  -1±) 


M 


cos  A  -  ■  ■  9  sin  A  cos  cp 

2ABZ 


+  U?  -  S  Sln  *  +  2k2  C0s  4  sin  *  +  N"  +  3  t  r2  '  <2-9 


and 


M 


'12 


N  z 


2AB 

1 


sin  A  cos  <J)  -  cos  A  sin  <J>  +  sin  4> 

zB 


M 


M 


2AB 


£•  2 


cB  R 


2AB 


(2.9.3c) 


where 


r  =  ± 

r  A 


(2.9.3d) 


The  metric  components  for  the  flat  are: 


G11 


1 


g22  =  N2  +  T jL?  (1  -  cos  A)2 


4AZB‘ 


(2.9.4a) 


3b) 


(2.9.4b) 
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and 

G12  =  0 

The  metric  components  for  the  middle  spiral  are: 
G11  =  0  > 


G22  =  N  +  2z 


n£  /±\ 


1_ 

2 


(2.9.4c) 

(2.9.5a) 

(2.9.5b) 


and 


G12  =  0 


(2.9.5c) 


The  Lagrangian  strain  tensor  components  are  defined  by 

eij  =  (Gij  '  gij)/2  1  (2.9.6) 

where  g..  are  the  metric  tensor  components  of  the  initial  flat 
ij 


plate : 


g .  .  =  6  .  . 

13 


(2.9.7) 


The  strain  tensor  components  for  the  edge  spiral  are: 


el  3  =  e2  3  “  e3  3  “  0 

N  z_  2 


el  1 

e2  2  = 
M 


c  R 

2 

M 


(2.9.8a) 

(2.9.8b) 


M  1 

2AB  I  (1  ’  COS  A)  ‘  4AB7  Sln  A  cos  ♦  +  SB? 

1 


4AB 


j  (1  -  cos  A)  sin  +  y  -  j  +  "B~  |-  r2  ,  (2.9.8c) 


and 


ei  2  - 


N 


M 


4AB 

1_ 

~2 


1  M 

sin  A  cos  <})  -  ^  +  —  (1 


-  cos  A)  sin  <f> 


2cB  R 


(2.9. 8d) 


Similarly,  for  the  flat, 

el  1  =  0 

e2  2 


Nz  _  1  1 

2  2  +  2 


M 


2AB 


(1  -  cos  A) 


(2.9.9a) 

(2.9.9b) 
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0 


(2.9.9c) 


and  for  the  middle  spiral, 


1_ 

2 


(2.9.10a) 


(2.9.10b) 


and 


e  1  2 


0 


(2.9.10c) 


It  is  readily  shown  that  due  to  deficiencies  of  the  mapping 


(such  as,  for  example,  that  discussed  in  section  2.6),  B  and  n  must 
be  large,  indeed,  much  larger  than  the  experimental  values  in  order 
that  the  middle  surface  strains  be  less  than  the  yield  point.  The 
experiments  of  Ricardo  (34)  have  shown  that  on  buckling,  the  change 
of  bending  energy  predominates.  These  experiments  confirm  Kennedy’s  (33) 
contention  that  the  deformation  is  largely  inextensional .  The  point 
of  view  is  therefore  taken  that  even  though  the  middle  plane  strain  of 
the  mapping  is  excessive,  the  primary  deformation  of  the  shell  is  due 
to  displacement  via  bending,  and  the  mapping  should  approximate  this. 

The  middle  plane  strain  energy  will  accordingly  be  computed  by  considering 
the  edge  and  middle  spirals  to  be  members  of  a  pinned  space  truss  (The 
geometry  of  the  space  truss  is  fixed  by  its  bending  deflections.), 
subjected  to  the  applied  axial  load.  Cox  (35)  has  noted  this  procedure 
as  a  possible  approximation  and  has  stated  that  the  flat  does  not 
participate  greatly  in  the  carrying  of  the  load.  Based  on  these  consider¬ 
ations,  the  membrane  terms  in  the  strain  tensors  will  be  dropped  from 
further  considerations.  The  parameter  N,  which  appears  in  the  bending 
terms  of  the  strain  tensor,  will  be  set  equal  to  one;  this  procedure  is 
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justified  if  the  axial  shortening  is  small  so  that  N  approximates 
unity. 

This  problem  with  the  middle  plane  strain  was  to  be  expected; 
one  has  but  to  note  Rayleigh's  comment  that  (35):  "We  can  bend  a 
piece  of  sheet  metal  easily  with  our  fingers  but  we  cannot  stretch 
it  noticeably".  Clearly,  the  large  deflection  mapping  would  have  to 
be  extremely  accurate  to  keep  the  strain  (which  must  be  compressive) , 
less  than,  say  0.0005  in. /in.  in  the  flat.  Indeed,  this  value  of  strain 
is  still  probably  much  larger  than  the  usual  post-buckling  values. 
Pogorelov's  (14)  procedure  approximated  a  surface  isometric  to  a  cylinder 
but,  as  in  this  thesis,  he  was  unable  to  maintain  geometric  continuity. 
Difficulties  with  the  middle  plane  strain  have  already  been  indicated 
in  the  discussion  of  Hoff,  et.  al. ,  (12,  1966). 


CHAPTER  3 


ENERGY  COMPUTATIONS  AND  THE  MINIMIZATION  PROCEDURE 


The  geometrical  considerations  of  Chapter  2  introduced  the 


parameters  n,  B,  k  and  N.  It  was  noted  that  c  was  subject  to  a 
geometric  constraint,  and  an  expression  c  =  c  (k,  n)  was  obtained. 
Thus,  if  c  is  prescribed,  and  n  is  chosen  as  an  independent  para¬ 
meter,  then  k  is  determined.  The  actual  procedure  shall  be  modified 
so  that  k  is  prescribed,  n  is  considered  as  independent,  and  hence 
c  is  determined.  A  consideration  of  the  magnitudes  and  continuity 
of  the  middle  plane  strains  indicated  that  while  the  mapping  may  be 
able  to  approximate  the  bending  deformation,  it  could  not  produce 
acceptable  middle  plane  deformations.  The  conclusion  was  reached  that 
the  middle  plane  deformation  would  have  to  be  accounted  for  in  an 
alternate,  more  reasonable,  fashion.  To  this  end,  N  was  set  equal  to 
one.  The  axial  shortening,  which  N  represented,  was  to  be  calculated 
as  a  function  of  the  remaining  parameters,  as  is  explained  below. 

The  principle  of  virtual  work  will  be  employed  to  determine 
B  and  n.  If  V  is  the  potential  energy,  then  it  is  required  that 


(3.1a) 


and 


(3.1b) 
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for  prescribed 
V  = 


k  and  c.  V  is  computed  from 
(U  -  W)/  -B°Lot 


where 


U 


and 


W 


U.  dS 
1 


(3.2) 


represents  the  area  upon  which  the  stress  vector  [i.e.,  (0,  0,  -a)] 
is  prescribed.  U  shall  be  calculated  as 


where  the  subscripts  denote  membrane  and  bending  energy,  respectively. 
The  regions  of  integration  are  shown  in  Fig.  24.  The  volume  of  this 
region  is : 


Further, 


and 


Volume 


BqLo  t 


2 


■  (?: 

j  |b2Ar|  t 

=  2BA2 

t3  (R/t)2 

(3.3) 

from  Fig.  24, 

d  =  L  cos  <f> 

9 

(3.4) 

Lq  =  (1  -  c) 

Lo  , 

(3.5) 

1 

1  =  B0  (B2  + 

0 . 25)  2 

(3.5a) 
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LO 


-x 


Fig.  24  Regions  of  Integration 


The  geometrical  shortening  e  (which  replaces  N)  is  obtained 
by  the  following  approximate  procedure.  Consider  a  section,  x  =  0, 

as  shown  in  Fig.  25.  The 


Fig.  25  Calculation  of  e 
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shortening,  e,  is  given  by 
e  =  1  -  cos  y 


where 


,-l 


MR 


=  sin  1  - —  (1  -  cos  A) 


Using  M  =  kA/sinA  and  Lq  =  BBg ,  Eq.  3.7  becomes 

Y  =  sln-l  K.U  -cos  A) 

2B  sin  A 


(3.6) 


(3.7) 


(3.8) 


The  expression  (see  Fig.  21) 

W  =  MR  (1  -  cos  A)  -  w  , 
was  taken  as  a  measure  of  the  maximum  lateral  deflection. 


(3.9) 


The  stresses  are  assumed  to  be  in  the  elastic  range  and 
are  given  by: 


and 


xx 


yy 


a 


xy 


E 


(1  -  v2) 


(e 

+ 

ve 

XX 

yy 

(e 

+ 

ve 

yy 

XX 

(1  +  v) 


xy 


3.1.  Calculation  of  the  Elastic  Energy  of  Bending  and  of  the  Potential 

of  the  Applied  Load 

The  bending  strain  tensor  components  for  the  edge  spiral 
are  (with  N  =  1)  : 


1  z  2 

:H  =  cRV 


(3.1.1a) 


1  Z.  2 

d  r  » 


4c  W-  R 


e22 


(3.1.1b) 
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and 


-1 


z 

—  r 


12  2cB  R 

Similarly,  for  the  middle  spiral, 

1 

^  /  z  \  d) 


'22 


./*  U  2 


Ug  is  calculated  from 


U 


B 


■  W 


v 


a  . .  e. .  dV 
ij  ij 


Expanding  Eq.  3.1.2  for  the  side  spiral,  there  results 
U. 


B 


=  4  f  (a  e  +a  e  +2a  e)dV 

l  J  xx  xx  yy  yy  xy  xy 


which,  after  substitution  from  the  stress-strain  laws,  becomes 


(3.1.1c) 


(3.1. Id) 


(3.1.2) 


U. 


B 


1  (  r  e  /  xx 

2  J  vLl  +  V  \  1  -  0 


+ 


r-EL  +  2e2 

1  -  v  xy 


2  E  v 


-r  -7-z  rT  e  e 

(1  -  v)z  XX  yyj 

An  exact  integration  over  the  region  shown  in  Fig.  24  was  excessively 
complicated;  since  the  method  is  approximate  in  nature,  the  regions  of 
integration  were  simplified  as  shown  in  Fig.  26. 


dV. 


(3.1.3) 


The  evaluation  of  Eq.  3.1.3  proceeds  as  follows: 


I  (  -X"  +  “ 

l  1  -  v  1 


+  2  eZ  dV  = 

-  v  xy  / 


1 

[  1  1 

fl  +  A.] 

1+  1  1 

^R? 

[  1-v  1 

[1  16B^  J 

2B7  J 

where 


y L'o 


x=x, 


//rdxdy  0  / 


x=x 


T  ,  x2 

fet2- 


dy 


4  2 

=  ~  cBA  (1 


2 

c)  R  , 


dx  dy 


'■ 
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and  also 


L 


e  e  dV 
xx  yy 


1 

4c2b"2 


P  J 


r^z  dV 


1  1  1  -  .2  3 

36  B2  A  C  4 


Fig.  26  Simplified  Region  of  Integration 


Combining  the  various  results,  the  elastic  bending  energy  of  two  edge 
spirals  is  given  by: 


U 


B 


9(1  -  P)  l  1  +  16B4  +  IP 


+ 


h ) 


1  -  c„„2  3 

- BA  t 


(3.1.4) 


The  elastic  bending  energy  for  the  middle  spiral  reduces  to 


U. 


B  2  ( 


V 


e  dV 

yy 


(3.1.5) 
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U 


B 


Et3 

24(1  -  v2 


y=d  x=+xn 


/  /  4>dxdy 

’ TRM  1  v=0  J  x=-x_ 


where 


<p  =  (-z)2/A2 


A  = 


2V 


,  and  z  =  y  -  d 


Upon  integration,  Eq.  3.1.5  reduces  to 


UT 


Ef 


Ar 

warz  (i  -  2c) 


'B  24(1  -  vO  R^y2  6 
which,  noting  Eqs.  2.7.1  and  3.4,  becomes 
UT 


Ef 


Y2  (1  -  2c) 


B 


18(1  -  v2)  coscj)  c 


However , 


COS<f>  = 


Ln 

\f(Ll  +  (B0/2)2) 

0 


B 


|/’(BZ  +  0.25) 


so  that  the  final  form  of  Eq.  3.1.5  is  1 

Et 3y 2  (B2  +  0 . 25) 2  (1  -  2c) 


U. 


B 


18(1  -  v2) 


B 


The  sum  of  the  bending  energies  therefore  yields 


U 


B 


EBA2t3  9j(l  -  v1) 

(B2  +  0.25)2  (1  -  2c) 

B2A2  ”  c 


i1  +W  +  IF)  1 


(3.1.6) 


-  c 


18(1  -  v2) 

(3.1.7) 


Equation  3.1.7  will  now  be  linearized  in  the  following  fashion.  From 


Eq.  3.8, 


y  = 


.  _i  k(l  -  cos  A) 

sm  1  — z  ~ - : - ; — 

2  B  sm  A 


•  -1  k(A2/2) 

=  Sin  2BA 

kA 

=  4B  • 


From  Eq.  2.5.6, 


c  = 


2A  -  2A 


A  -  M  sin  A 

- 


1  - 


5.2!  9.4! 


13 


(3.8a) 


+ 


r  _ 
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which,  using  Eqs.  2.6,9  and  2.6.10,  becomes 

_L  A  -  k(A/sinA)  sinA  5(1  -  k) 
°  A ^  _  2Ab~  '  2  A  7  A2 

5  9.4!  +  13.6! 


The  final  linearized  form  of  c  is 

(1  -  k)  n2 


(3.1.8) 


Equation  3.1.7  becomes 
U- 


B 

EBA2t3 


=  UB* 


1 


-f 


9(1  -  vz) 
1 


1  + 


I?B^  [  2B' 


)  ( 


18(1 


(1  -  k) n2 

_ _  (B2  +  0 . 25) 2  2  /  2  _  \ 

-  v2)  16  B4  \  (1  -  k)  n2  J 


-  1 


.  (3.1.9) 


The  potential  of  the  load  is 
-  (Bq/2  t  a)  L0  (1  -  cos  HO  = 
-2  (1  -  cos  HO  EEA2t3 


which , 


upon  normalization,  becomes 


W* 


/  oR  R  \  k2  7T 2 
\  Et  t  )  32  B2n2 


(3.1.10) 


3.2.  Approximation  to  the  Membrane  Energy 

If  the  spirals  are  treated  as  members  of  a  pinned  space 
truss,  then  a  typical  diamond  will  carry  an  axial  load  of  2AR  to  and 
a  typical  node  will  appear  as  in  Fig.  27. 
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Let  , 


C  =  compressive  (positive)  force  in  the  rib. 


and 


T  =  tensile  (positive)  force  in  the  middle  transition. 


The  equilibrium  equations  to  be  obeyed  are: 

2  •  (Y  component  of  C)  =  applied  load 

=  2AR  ta  ,  (3.2.1a) 


and 


2  •  (Radial  component  of  C)  =  radial  component  of  T. (3.2.1b) 
Some  necessary  geometrical  relations  are  obtained  from  Fig.  28: 


„  MR  (1  -  cos  A) 

a  =  Tan  1  - rrr - : - - - - 

MR  « i  n  A 


=  sin 


1 (1  -  cos  A) 


(3.2.2a) 


and 


Lq  cos  ¥ 

•  -l 
sm  1 


6 


1 


(3.2.2b) 
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The  Y  component  of 

C  is  C  sin  8  , 

and 

the  radial 

component 

is  (C  cos  8 )  sin  a.  The  radial 

component  of  T 

is 

T  sin  A. 

Equations 

3.2.1  therefore  become: 

2  C  sin  6  =  2  A  R  t  a 

and 


2  C  cos  8  sin  a  =  T  sin  A 
which  yield 

C  =  aAR  t/sin  8  , 


and 


T  =  2  C  cos  8  sin  a/sin  A 


(3.2.3a) 

(3.2.3b) 


Fig.  28  Approximate  Geometry  at  Node 
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Fig  0  29  Dimension  of  Members 


The  members  of  the  space  truss  are  approximated  as 


in  Fig.  29.  The  elastic  energy  of  the  side  spiral  is: 

C21 


U 


2AE 

1  1  1  q2ARt 

8  c  cos  <p  E  sin2  g 


A  =  2dt 


From  Eq .  3.2.2b, 


sin  g  = 


BBq  cos  y  . 


BB( 

~T 


=  COS  <j) 


Therefore , 


U 


1  1 


4  cBY  (0.25  4-  B2)2  ]  EBA2t3 


OR 


The  elastic  energy  of  the  middle  spiral  is: 

2 


u  =  i  /I)  11 

T  2  \  2  I  AE 


A  =  dt 


1  (0.25  +  B2)2 


32 


B 


k  /  oR  \ 
c  \  Et  ) 


sinA 


shown 


(3.2.4) 


c 


EBA2t3 


(3.2.5) 


. 
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Upon  normalization  with  respect  to  EBA2t3  and  linearization  of  c, 


Eqs.  3.2.4  and  3.2.5  become 
U 


2n2 


EBA^t 


27T 


1  (0.25  +  Bz) 


B 


(1 


_i _  /  or  y 

-  k)nz  \  Et  ) 


=  Uc*  ,  (3.2.6a) 


and 


U 


1_ 

2\  2 


EBAzt 


1_  (0.25  +  Bz) 
16  B4 


k 


(1  -  k)nz 


OR  \ 

Et  J 


=  UT*.(3.2.6b) 


3.3.  The  Minimization  Procedure 


Define : 


V*  =  U*  +  U*  +  U*  +  W* 
B  c  T 

The  governing  equations  become 
3V* 


3B 


=  0 


(3.2.7a) 


and 


9y* 

3n 


=  0 


(3.2.7b) 


for  prescribed  k  and  o.  In  expanded  form,  these  equations  are: 

1 

2 


3V*  _  ,  [ -1  1 

3g  k  l  /. n  5  b3 


{W  -h}{  (l--k)n^  - 

1 


,  (B2  +  0.25)  2  k2 

+  d  P  16 


2  „  \  (B2  +  0.25)  2k2/  -  2 

(1  -  k)n2  -  2  /  -  d  P  IiUl  -  k)n2  -  \ 

-  /OR  R  \  ^2  _k2_  /  (B2  +  0 . 25)  ~2  _  (B2  +  0. 25)  ^  k 

\  Et  t  j  8  B3n2  \  16BJ  4B3  j(l  -  k)n: 


aR 

Et 


I  (B2  +  0.25)  _  (B2  +  0 . 25) 2  _ 

+  2  |  ^3  p  j  (1  _  k)n2  ^Et 


OR  \2 


=  0 


(3.2.8a) 


and 


3V*  i  /  -i  .  1  ,1 

PT  =  b  1  +  I6B 4  +  IP 


-4 


(1  -  k) n3 


+  d 


(B2  +  0 . 25) 2 


B 


-4 


16  \  (1  -  k)n3j+  \  Et  t  j  8  B^n2  8 


aR  R  \  tt2  k2  1  (B2  +  0 . 25) 2 


B 


(1  -  k) 


*(f?) 
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(B2  +  0 . 25)  2 _ 1 


B 


(1  -  k)n- 


oR 

Et 


=  0 


(3.2.8b) 


where 


b  9(1  -  v2)  and  d  18(1  -  v2) 

The  parameter  k  should  be  very  close  to  one;  hence,  in  Eqs.  3.2.8, 

all  reference  to  k  other  than  in  the  (1  -  k)  form  will  be  dropped. 

These  equations  can  be  further  simplified  in  form  by  introducing  the 

following  symbols,  which  are  compatible  with  Fortran  coding:  \ 

1 


=  b 


A 

C  = 

E  = 

G  = 

H  = 

K  = 

M  =  - 


1 

4B- 


"  B*/! 


(B2  +  0.25)  2 


,  (B2  +  0 „ 25) 2 

-  d  VbS'  ~ 

(0.25  +  B2)"2 
16Bd 


D  = 
_1_ 

(0.25  +  B2)2 

4B3 


F  = 


2 

b 


(0.25  +  B2)  (0.25  +  B2)2 


BJ 

1  +  "~d  '4'  4-  — i- 
16B4  2B‘ 

TT2 


B 


OR  R 
Et  t 


8B 


7 


J  = 

L  = 


16B 

1 

f  oR  R  \ 

IT2 

,  Et  t  ) 

8B  6 

/  OR  \ 

2 

1  Et  ) 

J_ 

(B2  +  0. 

25)  2 

16B4 

1 

(0.25  + 

B2)2 

8B 


(0.25  +  B2)2 


B‘ 


Thus,  Eqs.  3.2.8  become 


9V* 

9B 

+  C 


=  A 


(1  -  k)  n 


(1  -  k)nz 
T  -  2 


-  1  1  +  BB  1  (1  -  k)n^ 


-  2 


n2(l  -  k) 


(  3.2.9a) 


and 


9n 


=  H 


-4 


(1  -  k) 


+  J 


-4 


(1  -  k) 


F(L  +  M)  _ 

+  K  +  (1  -  k)  ‘  ° 


(3.2.9b) 


It  was  previously  stated  that  the  equilibrium  equations  served  to 
define  n  and  B  for  prescribed  k  and  a  •  The  parameter  B  appears 
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in  Eqs .  3,2.9  in  a  very  complicated  fashion;  hence,  the  following 
procedure  was  adopted.  Solving  each  equation  for  (1  -  k)  there 
results 


2 (A  +  BB  +  C)  +  F(E  +  G) 

n2(A  +  2BB  +  2C)  -  D 


(3.2.10a) 


and 


(1  -  k) 


4H  +  4J  -  F(L  +  M) 

K 


(3.2.10b) 


Upon  equating  these  equations  one  obtains 

[2  (A  +  BB  +  C)  +  F(E  +  G)]K  \  _ 1 _ 

4H  +  4J  -  F(L  +  M)  °  J  A  +  2(BB  +  C) 

(3.2.11) 

Thus,  for  prescribed  B,  the  parameter  n  was  obtained  from  Eq.  3.2.11, 
and  k  was  determined  from  Eq.  3.2.10b.  A  program  was  written  to 
evaluate  these  equations  over  the  range  0.1  B  <_  0.75  using  B 
increments  of  0.01.  The  results  showed  that  the  neglect  of  k  in 
Eqs.  3.2,8  was  justified.  For  example,  when  aR/Et  =  0.1  and 

R/t  =  1000,  B  values  of  0.45,  0,40,  and  0.35  gave  values  for  (1  -  k) 

of  0.0045,  0.0048,  and  0.0052  respectively. 


CHAPTER  4 


RESULTS  AND  CONCLUSIONS 


4.1.  Results 

Because  of  the  failure  to  enforce  a  periodicity  constraint 
on  the  mapping  functions,  it  was  noted  that  the  parameter  k  would 
have  to  be  prescribed.  Rather  than  prescribe  k,  the  computations 
were  simplified  by  assigning  B  and  calculating  the  corresponding 
k  and  n.  Figures  30  to  38  contain  a  summary  of  the  results. 

The  parameter  B  appears  as  a  constant  in  the  various 
graphs.  Equation  3.2.9  was  not  valid  for  B  >_0.53  because  the 
solution  for  n2  became  negative.  The  wave  number  n  reached  a 
maximum  value  at  B  jl  0.31.  The  significance  of  these  points  is  not 
clearly  understood.  However,  it  is  noted  that  experiments  do  not, 
in  general,  show  B  >  0.5  if  the  buckling  occurs  periodically  over 
the  entire  surface;  the  model,  hopefully,  reflects  this  fact.  The 
exact  solution  for  the  periodic  buckle  pattern  would,  by  comparison  to 
experiment,  yield  B  <  0.5.  It  would  seem  that  the  procedure  used 
by  Hoff,  et.  al.,  (12)  (They  termed  their  method  of  solution  the 
von  Karman-Tsien-Legget  procedure)  produces  B  >  0.5  in  order  to 
prevent  the  axial  stress  due  to  deep  deflections  from  becoming  tensile. 
It  should  be  emphasized  that  the  axial  wave  length  of  Hoff,  et.  al. ,  is 
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approximately  six  times  that  indicated  by  experiment.  Recently, 

Jones  (16)  has  carried  through  the  von  Kdrman-Tsien-Legget  procedure 
for  prescribed  n;  it  would  be  of  interest  to  prescribe,  instead, 

B  jl  0.5,  vary  n,  and  to  then  take  note  of  whether  this  procedure 
would  force  a  larger  number  of  buckles  in  order  to  decrease  the  deep 
deflection  stresses.  Similarly,  it  might  be  of  value  to  revert  com¬ 
pletely  to  the  original  procedure  advanced  by  von  Karman  and  Tsien 
(8,  1941):  prescribing  both  n  and  B.  At  present  no  conclusion  can 
be  drawn  regarding  the  maximum  value  of  n;  however,  mathematically, 
one  expects  a  unique  point,  such  as  a  maximum  value  of  n,  to  occur 
on  a  aR/Et  =  constant  line;  but  where  it  should  occur  is  not  clear. 

It  was  observed  that  for  all  but  the  highest  R/t  considered, 
the  c  value  computed  using  B  =  0.5,  was  excessively  small  and 

did  not  correspond  to  experimental  values;  thus,  it  was  decided  that 
for  the  purposes  of  comparison  with  other  theoretical  and  experimental 
results,  the  results  for  B  =  0.31,  0.35,  0.40,  and  0.45  would  be 

considered . 

In  Eq.  3.2.8,  n2  and  t/R  were  not  related  so  that  all  refer¬ 
ence  to  n2  and  t/R  -  other  than  in  aR/Et  -  could  be  grouped  into 
the  dimensionless  parameter  n2t/R;  consequently  the  calculations 
were  carried  through  for  distinct  values  of  R/t.  This  is  in  contrast 
to  the  common  procedure  as  used  by  Hoff,  et.  al. ,  (12).  The  R/t 
values  considered  were  500  £  R/t  £  5000  in  increments  of  500. 
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Further  consideration  of  Eq.  3.2.9b  revealed  that  92V/8n2  =  0 

(computation  showed  32V/3B2  >0) ,  so  that  the  Hessian  determinant  is 
negative,  which  indicates,  unfortunately,  that  the  procedure  yields 
unstable  equilibrium  configurations.  It  is  possible  that  this  is  related 
to  the  essentially  unstable  results  of  Hoff,  et.  al.,  (12)  where  it 
was  noted  that  aR/Et  decreased  continuously  as  the  number  of  terms 
retained  in  the  w  function  increased.  The  nature  of  the  conclusions 
of  Hoff,  et.  al.,  will  be  further  discussed. 

The  comparison  of  computed  results  with  experimental  results 
is  made  difficult  by  the  nature  of  the  experimental  results  and  the 
nature  of  the  theoretical  assumptions.  By  way  of  example,  let  us 
consider  the  experiments  of  Thielemann  (21)  (see  Fig.  14).  Thielemann’s 
experiments  yielded  tier-type,  approximately  square,  buckles  and  there 
would  seem  to  be  little  doubt  that  the  unbuckled  material  served  to 
restrict  the  rotation  of  the  diamond  flats  and  so  decreased  the  axial 
shortening  from  that  which  one  would  expect  for  periodic  buckles  over 
the  entire  surface  as  assumed  in  theory.  Yet,  it  was  noted  when  discussing 
Fig.  14,  that  Almroth’s  and  Kempner's  axial  shortenings  compared  favor¬ 
ably  with  Thielemann's  experimental  results.  It  would  seem  that  this 
comparison  is  fortuitous.  Almroth  (11)  stated  that  his  n  was  less 
than  Thielemann’s  experimental  value;  indeed,  from  Hoff,  et.  al.,  (12), 

Case  3  at  e  =  e  yields  n  5.  This  low  value  of  n  should  serve 

cl 

t o  increase  e.  From  Hoff,  et.  al.,  it  is  further  seen  that  Almroth 
had  a  y  value  of  approximately  0.2,  that  is  to  say,  his  axial  wave 
length  is  about  five  times  larger  than  that  indicated  by  experiment . 
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Thus ,  a  wave  length  which  is  too  large,  has  combined  with  a  wave 
number  which  is  too  small,  to  yield  results  comparable  to  tier-type, 
approximately  square,  buckles,  whose  axial  shortening  has  been 
restricted  by  the  adjacent  material.  It  is  apparent  from  Yoshimura's 

(15)  approximate  formula  for  the  ideal  case 

2  2 
V  y 

6  = 


that,  for  y  =  1,  the  axial  shortening  is  severe.  Using  n  =  7 

(which  corresponds  to  Jones'  (16)  value  for  R/t  =  788  and 

oR/Et  =  0.1)  and  y  =  1,  there  results  e  =  0.0255,  which  bears 

no  comparison  to  the  e  value  (see  Fig.  14,  e  =e  ),  which  is 

cl  cl 

0.00078.  However,  replacing  y  by  0.2  reduces  e  to  0.00102  which 
is  comparable  to  the  classical  value.  As  previously  discussed,  this 
small  value  of  y  is  necessary  in  order  to  prevent  the  development 
of  axial  tension  stresses. 


Let  us  consider  the  series  of  graphs  in  Fig.  32a,  b,  c,  d, 
and  e.  These  graphs  show  that  oR/Et  decreases  as  n  decreases  and 
thus  indicate  the  unstable  nature  of  the  equilibrium  with  respect 
to  n.  It  is  noted  that  the  use  of  n  as  a  continuous  variable 
has  allowed  the  prescription  of  oR/Et  values  which  required  n  <  2 
in  order  to  obey  the  equilibrium  equations.  The  value  n  >_  2 
is  physically  necessary  in  order  that  the  solution  may  be  associated 
with  the  problem  of  a  closed  cylinder;  indeed,  stable  equilibrium 
would  seem  to  demand  n  >_  3,  since  n  =  2  is  a  hinge  and  would 

appear  to  be  naturally  unstable.  It  is  clear  that  n  increases  as 
R/t  increases  for  a  prescribed  oR/Et.  Further,  as  noted  above,  n 
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increases  with  oR/Et;  this  is  confirmed  by  the  experiments  of  Thielemann 
(21) .  Table  1  gives  a  comparison  of  the  results  of  this  thesis  with 
the  experiments  of  Thielemann  (21) .  For  a  given  n  ,  the  minimum  value 
of  the  experimental  load  obtained  by  Thielemann  is  quoted.  In  Table  1 
there  are  three  computed  values  of  n  given:  n2  is  the  non-integer 
value  associated  with  the  load  of  Thielemann,  while  n^  and  n^  are 
the  closest  integer  values.  For  oR/Et  =  0.1 

Table  1 

Comparison  to  Thielemann’s  Experiments 


R 

t 

Thielemann 

aR 

3 

aR 

n 

0.45 

0.40 

Et 

Et 

ni 

n2 

n3 

nl 

n2 

n  3 

788 

.119 

14 

.109 

.119 

.14 

7 

7.18 

8 

8 

8.8 

9 

.08 

6 

8 

.109 

13 

.109 

7 

8.4 

.14 

8 

9 

.08 

6 

7 

.100 

12 

.100 

6.55 

7.8 

.  109 

7 

8 

.06 

5 

6 

.0763 

11 

.0763 

5.96 

7 

.08 

6 

8 

.056 

5 

5 

.066 

10 

.066 

5.35 

6.5 

.08 

6 

7 

' 
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Jones  (16,  1966)  obtained  n  =  7.  Mo  comparison  can  be  made  with  the 
results  of  Hoff,  et.  al.,  (12)  since  their  solution  lies  below  the 
values  of  the  loads  quoted.  It  is  observed  that  the  experimental  values 
of  n  exceed  the  calculated  values.  Further  comparison  is  possible 
with  the  very  refined  experiments  of  Almroth,  Holmes,  and  Brush  (31). 

For  an  R/t  value  of  850,  these  investigators  obtained  post-buckling 
loads  of  0.08  to  0.1,  with  the  number  of  buckles  varying  from  8  to  12. 
For  aR/Et  =  0.08,  the  method  of  calculation  used  in  this  thesis 

yields  n  =  8.15,  8.1,  7.5  and  6.3  corresponding  to  B  values  of 

0.31,  0.35,  0.40  and  0.45.  Similarly,  for  aR/Et  =  0.1  the  values 

obtained  were  n  =  9.15,  9.1,  8.4,  and  7.  Hoff,  et.  al.,  (12) 

obtained  a  minimum  aR/Et  of  0.0427.  The  value  of  y  was  about  0.17 
(this  is  to  be  compared  to  the  experimental  value  of  1) ,  and  for 
R/t  =  1000,  the  value  of  n  was  close  to  three.  The  calculations  of 
this  thesis,  applied  to  the  same  case,  yielded  n  =  6.6,  6.5,  5.9, 
and  4.9  for  B  values  of  0.31,  0.35,  0.40,  and  0.45. 

The  fact  that  n  decreases  to  a  value  less  than  two  would 
seem  to  correspond  to  the  results  of  Hoff,  et.  al. ,  (12).  Suppose  one 
applies  the  results  of  Hoff,  et.  al. ,  to  a  specific  R/t;  then  it  is 

required  that,  in  order  for  the  solution  to  correspond  to  that  for  a 

closed  cylinder,  there  exist  the  relation 

p  =  n^t/R  >_  2^t/R  . 

Thus,  if  the  calculated  value  of  n  is  less  than  this  bound,  the 
equations  are  giving  a  value  of  n  less  than  two  and  the  solution  can 

no  longer  be  identified  with  the  problem  of  a  closed  cylinder.  In 
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the  light  of  the  fact  that  the  computations  of  Hoff,  et.  al. ,  give,  in 
this  authors  opinion,  decreasing  values  of  n  as  the  number  of 
parameters  is  increased,  let  us  return  to  the  paper  of  Hoff,  et.  al. 

The  synopsis  states  in  part:  "The  results  of  the  computations  as  well 
as  theoretical  considerations  indicate  that  the  solution  for  long  shells 
of  the  von  K^rman-Donnell  equations  with  the  aid  of  the  von  Karman- 
Tsien-Legget  procedure  leads  in  the  limit  to  the  trivial  solution  in 
which  the  amplitudes  of  the  displacements  tends  to  zero,  the  number 
of  waves  around  the  circumference  tends  to  infinity,  and  the  average 
axial  compressive  stress  capable  of  maintaining  equilibrium  in  the 
post-buckling  state  tends  to  zero."  This  would  not  appear  to  agree 
with  the  ideas  presented  above,  and  in  this  authors  opinion,  cannot 
be  deduced  from  the  von  Kdrm£n-Tsien-Legget  approach.  It  is  the  opinion 
of  this  author  that  Hoff,  et.  al. ,  used  the  following  line  of  reasoning. 
They  first  required  that  the  values  of  the  strains  computed  from  the 
Donnell  equations  approach  the  exact;  this  required  that  n  tend  to 
infinity.  That  n  tended  to  infinity  then  agreed  with  the  result 
that 

y  =  1  /I  =  2ttR/  2nl 

y  x  x 

decreases  to  zero.  Next,  since  n  has  been  allowed  to  tend  to  infinity 
and  since  the  von  Karmdn-Tsien-Legget  procedure  yielded  that 

ri  =  n2t/R 

approach  zero,  then  it  was  required  that  t/R  go  to  zero.  It  is  noted 
that  here  Hoff,  et.  al.  ,  used  the  argument  n  >_  2;  however,  n  tending 
to  infinity  is  needed  for  the  exactness  of  the  geometric  relations. 

Last,  it  was  observed  that  the  von  Karmdn-Tsien-Legget  procedure  gave 


* 
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aR/Et  decreasing  to  zero.  Thus  the  von  K^rmdn-Tsien-Legget  approach 
does  not  yield  the  solution  that  n  goes  to  infinity;  rather,  this 
would  appear  to  be  an  arbitrarily  introduced  constraint  to  force  exactness 
in  Donnell's  geometric  relations.  It  should  be  observed  that  allowing 
1  to  increase  and  n  to  decrease  can  give  the  solutions  of  p  and  n 

X 

tending  to  zero;  this  argument  is  independent  of  t/R  and  uses  only 
the  results  of  the  von  Kdrmdn-Tsien-Legget  procedure.  Certainly,  Fig.  15a 
does  not  indicate  that,  as  oR/Et  decreases  to  zero,  eR/ t  also 
decreases  to  zero,  as  would  be  required  if  n  increased  to  infinity. 

Let  us  further  consider  the  conclusion  of  Hoff,  et.  al.,  that:  "In 
the  limiting  case  n  approaches  infinity,  the  amplitude  of  the  dis¬ 
placements  tends  to  zero,  the  von  K^rmdn-Donnell  equations  are  rigorously 
valid,  the  t/R  approaches  zero,  and  the  value  of  average  axial  stress 
capable  of  maintaining  equilibrium  in  the  post-buckling  state  is  zero." 
Since,  in  their  opinion,  n  tends  to  infinity,  it  would  seem  that  the 
authors  are  unable  to  specify  t/R;  yet,  the  von  K3rman-Donnell 
equations  are  written  for  a  specific  t/R.  The  conclusion  would  there¬ 
fore  seem  to  be  that  t  =  0.  With  t  =  0,  the  governing  partial 
differential  equations  are  no  longer  the  von  K5rmdn-Donnell  equations 
subject  to  we  C8  but  could  be,  say, 

Got6  "  ga6  =  ° 

the  solution  of  which  allows  discontinuities  of  slope  and  curvature 
and  whose  solution  does  not  demand  that  n  go  to  infinity.  The 
ideal  diamond  shape  is  a  valid  solution  to  the  new  system  of  equations 


with  n  >  2. 


. 
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The  simple  model  proposed  in  this  thesis  has  aided  in  explaining 
the  essentially  unstable  results  of  Hoff,  et.  al.,  (12).  Previously, 
it  was  noted  that  32V/9B2  was  positive;  thus,  if  rather  than  employing 
the  equation  3V/3n  =  0,  one  instead  followed  the  method  of  Madsen 
and  Hoff  (27,  1965)  and  prescribed  n  (at  the  values  found  in  this 
thesis),  then  stable  configurations  would  be  obtained. 

Figures  33a  and  33b  give  the  variation  of  axial  shortening 
for  B  =  0.45  and  0.40.  The  value  of  e  is  calculated  from  a 

linearized  form  of  Eq.  3.6  which  corresponds  to  Yoshimura's  (15) 
approximation  for  the  perfect  diamond  shape  and  therefore  probably 
gives  an  e  value  which  is  too  large.  However,  purely  geometric 
considerations  applied  to  the  ideal  experimental  periodic  buckle 
pattern,  which  is  extremely  flat,  would  reasonably  indicate  shortenings 
close  to  the  values  as  calculated  above.  As  previously  mentioned,  the 
ideal  seldom  happens  and  so  it  is  difficult  to  compare  theoretical 
and  experimental  values  for  axial  shortening.  Hoff,  et.  al.,  (12) 
obtained  values  comparable  to  the  ideal  diamond  shape  for  a  given  y 
and  n;  but,  these  results  cannot,  reasonably,  be  compared  to  experiment 
since  their  axial  wave  length  was  about  six  times  the  experimental 
value . 


Figure  34  gives  an  indication  of  the  maximum  lateral  deflection 
for  the  cases  of  B  =  0.45  and  0.40.  It  is  clear  that  the  deflections 
and  rotations  are  extremely  large  and  as  such,  careful  consideration 
must  be  given  to  including  the  non-linear  effects  of  u  and  v  in  the 
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strain  tensor  components.  With  regard  to  the  large  rotations,  Hoff, 
et.  al.,  stated:  "Indeed,  use  of  the  stress  function  insures  equili¬ 
brium  in  the  tangent  plane  only  if  the  rotations  are  small."  This 
does  not  agree  with  Langhaar  (35)  who  states  the  stress  function  can 
be  allowed  for  almost  inextensional  buckling.  There  would  appear  to 
be  a  matter  of  interpretation  here.  It  is  assumed  that  they  mean  the 
tangent  plane  to  the  original  cylinder.  If,  indeed,  this  is  so,  then 
Hoff,  et.  al. ,  have  failed  to  recognize  that  x  and  y  were  used  as 
Lagrangian  co-ordinates  (otherwise,  the  volume  integrations  would  have 
been  excessively  difficult) ,  and  the  equilibrium  equations  are  therefore 
written  with  respect  to  a  local  reference  frame  based  upon  the  deformed 
tangent  plane,  and  not  the  original. 

The  dependence  of  c  upon  aR/Et  is  given  in  Fig.  35  for 
B  =  0.50,  0.45,  0.40,  0.35,  and  0.31;  there  is  no  R/t  effect  in  c. 

As  indicated  at  the  beginning  of  this  chapter,  B  =  0.50  yields  a 

value  for  c  which  is  too  small  and  which  would  keep  all  but  the 
highest  R/t  values  plastic.  This  is  undesirable,  since  a  largely 
elastic  configuration  was  assumed.  Further,  if  one  notes  that  the 
total  spiral  width  of  one  buckle  is 

2L  =  2(2cRA)  =  2c(2RA)  , 

and  that  the  maximum  buckle  width  is  2RA,  then  it  is  necessary  that 
2c  <  1 

and  a  more  reasonable  limit  would  seem  to  be 
2c  <  0.50 

Based  upon  these  considerations,  B  =  0.35  and  0.31  give  c  values 

which  render  the  geometry  extremely  approximate  in  nature. 
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Pogorelov  (14)  gave  the  maximum  bending  strain  as  3t/R 
(virgin  state  of  the  material  was  the  undeformed  cylinder) ;  he  obtained 
a  minimum  post-buckling  load  of  oR/Et  =  0.15  with  B  <_  0.455.  For 

purposes  of  comparison  with  Pogorelov' s  results,  Fig.  30  gives  the 
maximum  normal  strain  for  oR/Et  =  0.15  and  0.10.  It  is  clear 

that  the  maximum  bending  strain  calculated  in  this  thesis  exceeds  that 
of  Pogorelov.  The  discrepancy  may  be  partly  explained  by  the  fact  that 
Pogorelov  used  the  classical  wave  number.  Let  us  consider  the  case 
of  oR/Et  =  0.15  and  R/t  =  1000.  Pogorelov  computes  n  from 

n( t/R) 2  =  0.9  , 

which  yields 

n  =  28.5  ; 

while  the  results  of  the  present  thesis  give  for  B  =  0.4  the  value 

n  =  11.3 


e 

max 


R/t 


.30  .35  .40  .45  .50  B 

Fig.  30  Maximum  Bending  Strains 
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It  is  therefore  clear  that  the  lower  n  value  of  this  thesis  yields 
higher  bending  strains.  Figure  30  indicates  a  minimum  value  of  strain 
at  B  ^  0.36.  Figure  30  further  indicates  that  B  is  a  very  important 
parameter;  a  seemingly  small  change  in  B  results  in  a  considerable 
change  of  the  maximum  normal  strain. 

Tests  were  performed  using  cylinders  of  R/t  =  950  (the 

proportional  limit  of  the  material  was  approximately  68000  p.s.i.), 
and  plastic  strain  was  noted  on  the  ribs  as  well  as,  of  course,  at 
the  singular  points.  While  discussing  the  topic  of  plastic  strains, 
it  is  useful  to  consider  Figs.  36a  and  36b  which  give  the  elastic  energy 
density  for  B  =  0.45  and  0.40  (E  was  taken  as  30. 106  p.s.i.).  For 

yield  point  values  of  30,000  p.s.i.  and  60,000  p.s.i.,  the  limiting 
values  of  elastic  energy  density  are  15  p.s.i.  and  60  p.s.i.,  respectively. 
Table  2  summarizes  the  minimum  value  of  R/t  as  needed  to  maintain 
an  average  fully  elastic  energy  density  at  aR/Et  =  0.10. 


Table  2 

Critical  R/t  values 


Yield 

Point 

B 

0.45 

0.40 

R/t 

30000 

2150 

1900 

60000 

1100 

950 

For  high  R/t,  the  strains  are  largely  elastic.  Figure  37,  applied  to 
a  case  of  high  R/t  (this  insures  that  the  correct  stress  strain  law 
has  been  used)  ,  indicates  that  at  the  approximate  experimental  value 


■ 


' 
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of  aR/Et  =  0.1,  the  elastic  energy  is  97%  due  to  bending.  This 

would  imply  that  Yoshimura's  (15)  comment  stating  that  the  bending 
and  membrane  energies  are  comparable  is  not  accurate.  This  would 
further  justify  Pogorelov's  (14)  procedure  and  would  tend  to  verify 
Kennedy's  (33)  suggestion  that  the  post-buckling  configuration  be 
treated  as  largely  inextensional .  Indeed,  computations  which  neglected 
the  membrane  energy  were  performed  for  R/t  =  500,  and  it  was  found 

that  there  was  no  significant  change  in  the  results.  Figure  36b 
contains  plots  of  the  pre  -buckling  energy;  these  plots  indicate  that 
the  elastic  energy  increases  greatly  on  buckling.  Von  K^rmdn  and 
Tsien  (8)  had  indicated  that  a  release  of  elastic  energy  on  buckling 
served  to  explain  the  rapidity  of  the  buckling;  von  Kerman's  idea  may, 
however,  apply  to  local  buckling.  It  must  of  course  be  remembered 
when  commenting  upon  Yoshimura's,  and  von  Kerman's  statements,  that 
the  basis  for  the  criticism  is  an  acknowledged  highly  approximate 
model . 


It  was  desired  to  test  what  part  the  flat  played  in  carrying 
the  axial  load;  it  will  be  recalled  that  the  load  was  assumed  to  be 
carried  by  a  pinned  space  truss  which  had  the  edge  and  middle  spirals 
as  members.  Figures  38a  and  38b  give  plots  of  the  rib  compressive 
stress  for  B  =  0.45  and  0.40.  For  low  R/t  values  the  compressive 
stress  seems  excessive  and  this  would  indicate  that  the  flat  does 
participate  in  carrying  the  load.  If  one  could  devise  a  loading  system 
that  could  force  an  elastic  periodic  buckle  pattern,  it  would  be  of 
interest  to  employ  back-to-back  strain  rosettes  to  determine  how  the 


!  . 
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load  is  carried  through  a  typical  buckle.  An  experimental  study 
(using  steel  shells  of  R/t  =  950  and  proportional  limit  of  68000  p.s.i.) 

was  initiated  with  this  aim,  but  it  was  necessary  to  abandon  the 
experiments  since  it  was  impossible  to  obtain  an  elastic  configuration. 

The  experiments  used  a  buckle  width  of  n  =  8  as  determined  from 
random  buckle  patterns.  It  is  therefore  clear  that  if  one  uses  steel 
shells  in  such  experiments,  then  higher  R/t  ratios  are  needed. 

A  particularly  interesting  feature  of  the  model  presented  in 
this  thesis  is  the  nature  of  the  post-buckling  V  surface  envelope, 
as  determined  by  prescribing  aR/Et  and  B  and  solving  for  n  and 
k  in  the  manner  previously  described.  The  envelope  is  a  surface  of 
translation  and  is  generated  by  translating  the  line 


-1 

[  /  i  +1+1 

\  .  (B2  +  0 . 25) 2  1 

E 

9(1  -  v2) 

\  1  +  16B4  +  2B7 

/  +  16B4 

2 (R/t) 2 

parallel  to  the  n  axis  in  B  -  n  -  V  space.  Figure  31  gives  a 
diagram  of  this  surface.  Also  shown  are  the  oR/Et  =  constant  lines 
projected  onto  the  n  -  B  plane.  The  neutral  nature  of  the  V  envelope 
with  respect  to  n  is  illustrated  by  Fig.  31;  the  diagram  shows  that 
there  exist  equilibrium  states  characterized  by  constant  B  and  V 
values  and  in  which  n  decreases  as  oR/Et  decreases.  This  conservation 
of  potential  is  not  to  be  confused  with  Tsien’s  [(36,  1942)  and  (37, 

1947)]  proposed  buckling  criterion.  Tsien  considered  V  conserved 
between  the  pre-buckling  and  post-buckling  configurations  for  the 
separate  cases  of  aR/Et  prescribed,  or  e  prescribed.  Figure  31  is 
strictly  for  the  post-buckling  configuration.  Fung  and  Sechler  (38) 
have  stated  that  Tsien’s  criterion  is  not  based  on  sound  principles  and 
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have  given  an  example  where  it  fails  completely.  However,  in  the 
present  case,  V  in  the  post-buckling  state  is  conserved  for  constant 
B  as  an  analytical  result,  as  opposed  to  an  enforced  constraint. 


V 


Fig.  31  B-n-V  Space,  R/t  =  constant 


It  is  noted  that  in  this  thesis  no  distinction  has  been  made  between 
the  equilibrium  conditions  for  aR/Et  prescribed  and  e  prescribed. 
This  is  not  necessary,  since  one  expects  a  unique  functional  relation¬ 
ship  between  e  and  a,  say,  e  =  e(o) ;  and  further  one  expects 
the  functional  relationship  can  be  inverted  as  would  be  done  by 
prescribing  e  and  determining  a  =  o(e). 

4.2.  Conclusions 

The  author  feels  that  the  model  used  in  this  thesis  has 
served  to  illustrate  the  following  points: 
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1)  The  model  has  shown  the  inherent  instability  of  the 
system  if  the  only  restriction  on  n  is  3V/3n  =  0;  thus,  some 

additional  criterion  is  necessary  to  determine  the  value  of  n.  As 
previously  mentioned  the  lower  bound  to  n  must  be  at  least  three. 
By  physical  argument,  the  prescription  of  a  aR/Et  value  which 
necessitates  n  =  2  can  be  ruled  out,  since  this  load  cannot  be 

held  in  a  stable  equilibrium  configuration,  and  if  n  <  2,  then  the 
solution  no  longer  corresponds  to  a  closed  shell  problem. 


The  determination  of  a  criterion  to  fix  the  value  of  n 


is  of  the  upmost  importance.  The  procedure  used  by  Jones  (16,  1966) 
and  by  Madsen  and  Hoff  (27,  1965)  is  considered  to  be  unacceptable. 
They  used  the  von  Karman-Tsien  model;  the  form  of  the  w  deflection 


function  was 


w  = 


i=0  1,1 


2i7TX 

lx  ly  “2i,0  lx 


„  A  irx  iry  ,  _  . 

E  A.  .  cos  - -  cos  - — ^  +  E  A0 . 


cos 


The  unknowns  A.  .  ,  n,  and  y  were  determined  from  the  equations 

1  J  J 


3  V 


3A.  . 

i.i 


3  V 


3A0 .  n 

2i  ,0 


=  0 


=  0 


i  >  1 


F  -  ° 

3y 


and 

n  =  prescribed  =  ng 

The  value  of  Ag i  was  determined  from  the  periodicity  constraint  on  v. 


It  is  most  unsatisfying  to  prescribe  n  without  making  it 
subject  to  some  constraint;  further,  if  one  wishes  to  prescribe  n  then 
it  also  seems  reasonable  to  prescribe  y  since  experiments  indicate 
that  y  is  also  largely  constant  in  the  post-buckling  era.  It  is 
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therefore  suggested  that  the  procedure  of  Madsen  and  Hoff  (27)  be 
modified  to  include  the  prescription  of  both  n  and  y  but  that  there 
be  imposed  the  constraint  conditions 


3 V 
3n 


n=n0 


=  0 


and 


3V 

3y 


y=yo 


=  0 


The  constraint  conditions  may  be  considered  as  equations  for,  say, 

A22  and  A3 3.  The  periodicity  constraint  on  v  may  be  used  to  determine 
A} 1 .  The  final  set  of  equations  to  be  solved  would  be 
3V 


3A. 


=  0 


i  >  3 


and 


3V 


8A2i,0 


=  0 


The  results  of  such  computations,  when  compared  to  the  results  of 
Madsen  and  Hoff  (27) ,  will  indicate  the  effect  of  the  additional 
constraints  on  the  solution. 

2)  The  use  of  a  smaller  number  of  variables  has  eased  the 
analysis  and  has  yielded  results  comparable  to  other  current  methods. 

The  parameter  B  has  been  shown  to  be  most  important.  It  is  of 
interest  to  note  that  the  approximation  of  middle  plane  strain  energy 
used  in  this  study  indicates  values  of  B  much  closer  to  experimental 
results  than  are  the  values  of  Hoff,  et.  al.,  (12).  It  has  been  noted 
that  the  large  axial  wave  lengths  of  Hoff,  et.  al.,  are  necessary  to 
prevent  the  development  of  axial  tension  stresses  due  to  deep  deflections. 

3)  The  results  indicate  that  the  middle  surface  deformation 


is  largely  inextensional .  Further,  it  has  been  shown  that  under  the 
assumption  of  periodic  buckling,  the  post-buckling  elastic  energy  exceeds 
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the  pre-buckling  values  at  the  same  aR/Et.  Indeed,  the  results 
indicate  that  for  aR/Et  =  0.1,  97%  of  the  elastic  energy  is  due 

to  bending. 

4)  The  computations  indicate  that  if  periodic  buckling 
is  assumed  to  occur  over  the  entire  surface,  then  plasticity  effects 
must  be  considered.  For  example,  in  Fig.  30,  the  case  aR/Et  =  0.1 

requires,  if  the  yield  point  is  60,000  p.s.i.,  that  the  R/t  ratio 
should  exceed  3000  to  ensure  an  essentially  elastic  configuration. 
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Fig.  32 .  a 


UR/Et  versus  n 


B  =  0.  50 
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Fig.  32. b 


ITR/Et  versus  n 


B  =  0.45 
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Fig.  3  2 .  c 


tfR/Et  versus  n 


B  =  0. 40 
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NUMBER  OF  BUCKLES  -  n 

Fig.  32. d  -  \JR/Et  versus  n  -  B  =  0.35 
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Fig.  32. e 


tfR/Et  versus  n 


B  =  0.  3  I 
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AXIAL  SHORTENING  -  e 

Fig.  33.a  -  TR/Et  versus  e  -  B  =  0.45 
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Fig.  33.  b 


tfR/Et  versus  e 


B  =  0.40 
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MEASURE  OF  MAXIMUM  DEFLECTION  >  UNITS  -  t 


Fig.  34 


TR/Et  versus  MAXIMUM  LATERAL  DEFLECTION 
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Fig.  35 


U'R/Et  VERSUS  c 
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Fig  .  36. a 


TR/Et  VERSUS  ELASTIC  ENERGY  PER  UNIT  VOLUME 


B=  0 . 45 
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Fig.  36. b 


g R/Et  versus  ELASTIC  ENERGY  PER  UNIT  VOLUME 


B= 0.40 
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PERCENTAGE  BENDING  ENERGY 


Fig.  37 


U’R/Et  versus  PERCENTAGE  BENDING  ENERGY 


-108- 


Fig.  38. a 


tfR/Et  VERSUS  RIB  AXIAL  COMPRESSIVE  STRESS 


B  =  0 . 4  5 
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Fig.  38. b 


tfR/Et  VERSUS  RIB  AXIAL  COMPRESSIVE  STRESS 


B=  0.40 
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